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Chapter  I 
IM'RODUCTIOK 

Inundation  of  river  flood  plains  by  occasional  flood 
discharges  has  IjQeii  a  threat  to  the  life  and  property  of  L:an 
since  the  daim  of  civilization.    Man  has  long  sought  means " 
of  predicting  floods  and  averting  their  danger.    Modern  man's 
need  for  quantitative  IrnTorination  concerning  flood  flo\?s 
becomes  evc;n  more  pref^siDg  as  he  alters  r*atural  ijatcrways 
^^ith  Btructin^as^  realigmients,  and  diveroionSc 

A  major  problem  continually  faced  by  highway  designsrB,  urban 
plamerS;,  and  ^^atarshea  KanagcLieTifc  englDesrs  is  that  of  determining  the 
frequency  of  psak  flood  discharges  from  stream  and  river  basins e  For 
large  river  basins  (in  exC'-^ss  of  100  square  miles)  flood  frequency  in- 
formation is  gensrally  available  by  virtue  of  long  term  records  of  pre- 
cipitation_j  stream  f lo'? ,  etc,  vhich  have  been  collected  by  various 
public  agencies  (including  the  U.  Sc  Env i r onsie nt a  1  Science  Services 
Administration,  Geological  Stn'vey,  Bureau  of  Reclamtion,  Agriculture 
Research  Service,  and  others )«,    For  watersheds  smaller  than  100  square 
miles  in  area  there  is  generally  a  shortage  of  hydrologic  datao  Cnly 
within  the  past  feiv  years  has  significant  research  been  directed  toward 
a  study  of  flood  frequencies  on  basins  of  this  size.  ■  '  . 

The  "Drainage  Correlation  Research  Project ",  which  was  initiated  by 
the  Department  of  Civil  Engineerirg  and  EngineeriDg  M^^chanics,  Montarja 
State  University,  Bozeinan,  Montana,  in  I963,  addresses  itself  to  the 
problem  of  predictirg  the  frequency  with  which  a  peal:  discharge  of  given 
magnitude  luay  be  expactel  on  small  watersheds  in  Montana c    This  investi- 
gation is  sponsored  by  the  Montana  State  Highway  Department  and  the  U,  S, 
•Bureau  of  Public  Roads ^    Kie  work  is  being  done  under  the  direction  of 


Theodore  T.  Williaras,  Associate  Professor  of  Civil  Engineering  and  En- 
gineering ^^2chanics  at  Montana  State  University. 

The  Drainage  Correlation  Research  Project  is  a  two  phase  study: 
a)  to  determine  flood  frequency  relationships  on  watersheds  smaller  than 
100  square  miles  from  long-term  climatological  data  already  in  existence; 
and  b)  to  make  comprehensive  hydrologic  studies  of  four  small  v;atersheds 
in  eastern  Montana.    The  watersheds  selected  for  study  are  Bacon  Creek 
in  Wheatland  County,  Duck  Creek  in  Prairie  and  McCone  Counties,  Hump 
Creek  in  Sweet  Grass  County,  and  Lone  hhn  Coulee  in  Pondera  County,  Con- 
tinuous records  of  climatic  factors  and  streamflow  are  being  obtained  from 
these  four  basins.    A  variety  of  peak-flow  prediction  techniques,  which 
utilize  data  such  as  that  being  collected,  are  being  examined. 

Among  the  techniques  being  considered  under  the  second  phase  of  the 
Drainage  Correlation  Research  Project  is  one  which  was  developed  by  the 
Soil  Conservation  Service,  U.  S,  Department  of  Agriculture  (SCS),  The 
SCS  method  possesses  n^ny  desirable  features,  and  seems  to  have  consid- 
erable potential  as  a  tool  for  predicting  flood  frequencies.    The  study 
reported  in  this  report  analyzes  the  method  developed  by  the  SCS;  describes 
the  use  of  a  related  computer  program;  tests  the  program's  ability  to  re- 
produce actual  hydrographs  from  a  selected  watershed;  and  evaluates  the 
possible  utility  of  this  method  on  I-fontana  watersheds. 

The  investigation  reported  herein  consisted  of  the  analysis,  using 
the  SCS  method,  of  a  single  rain-caused  runoff  event  which  occurred  on 
Duck  Creek,  the  largest  of  the  four  watersheds    being  studied  by  the 


DrainsgG  CorrcO-atir-n  Research:  Project- «,    Riek  Crfcol:  is  an  aphair-jral  s'brefiisi 
typical  rf  most  caGtern  Kontam-  streams  iJhicb  dr£ii>  srall  uatershadSe  It 
i£?  dry  ?no3t  of  the  t;lne,  but  occasionally  has  aiiri"ac:u  f lo-.f  after  a.  rain- 
caused  runoff  event  or  durix:g  spring  •snox-zmslte    In  the  four  years  for 
which  data  have  been  obtained  at  Tack  Creek,  there  have  been  only  a  few 
runoff  event B  of  consequence and  only  an  evept  uliich  occurred  on  June  l6p 
1965  lends  itself  to  an  analysis  by  the  SOS  method c 

A  number  of  events  would  need  to  be  analyzed  before  conclusive  state- 
laentQ  as  to  the  validity  of  the  SOB  method  could  ba  madeo  Mevertheless_, 
it  is  believed  that  the  results  of  the  analysis  of  this  one  storra  will  be 
a  valuable  contribution  in  estimtiog  the  applicability  of  the  method  to 
Montana  watersheds c  ;  •• 


Chapter  II 

LITERATURE  REVIEW 

Scientific  hydrology  is  a  relatively  new  area  of  study.  Appli- 
cation of  the  scientific  approach  to  hydrologic  problems  has  had  its 
greatest  growth  in  this  century,  and  especially  in  the  last  thirty 
years. 

Hydrologic  oieasuretnents  were  not  considered  to  be  of  much  import- 
ance during  the  early  history  of  the  United  States.     Population  centers 
developed  along  the  rivers  and  lakes  and  so  no  shortage  of  water  was 
experienced  until  the  westward  expansion  into  the  arid  and  semi- arid 
regions.     The  beginning  of  systematic  collection  of  hydrologic  data  in 
this  country  can  probably  be  taken  to  be  in  I888,  when  the  U.  S.  Geo- 
logical Survey  under  the  direction  of  Frederick  H.  Newell  set  up  its 
first  river-measurement  station  on  the  Rio  Grande  R^ver  at  Embudo,  New 
Mexico. 

Modern  Hydrology 

Since  about  1930?  the  increase  in  the  amount  of  printed  information 
made  available  in  the  field  of  hydrology  has  increased  very  rapidly. 
This  indicates  something  of  the  increased  need  there  has  been  in  recent 
years  to  obtain  more  data  and  develop  better  prediction  methods. 

The  study  of  hydrology  can  be  subdivided  into  the  two  general  areas 
of  stochastic  and  deterministic  hydrology.     Deterministic  hydrology  can 
be  further  separated  into  physical  or  analytical,  dynamic  and  parametric. 

Stochastic  hydrology  attempts  to  predict  future  events  or  reconstruct 
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past  events  based  on  tha  statistical  properties  of  the  kno",;n  record c 
White  (1967)  states  that  if  the  value  that  e.  variab2.e  has  iHiXJlies  an 
eleraent  of  chance^  then  it  is  a  stochastic  variable  ^    Thus  if  ths  record 
avallabla  for  so?i9  variable  (the  value  of  uhinh  is  ran'Ioa  in  nature  or 
stochastic)  is  restricted  in  tina,  then  the  predictecl  future  projection 
of  thie  record  is  hooa  by  statistical  Ksans-e 

Physical  {ai:;ajjtical )  hydrcdo^y  addresses  itself  to  particular 
specialized  problems  of  ths  hydrolcgic  procGGs ,    It  doss  not  try  to  supply 
any  ariBuers  to  paak  discharge,  aiirmal  yiGld_;  or  frequency  studies «  In 
a  sense  this  branch  of  hydrology  iB  i;ure  research  in  that  it  seeks  to 
establish  relat ionahips  betvjeen  variables  operatiri^  in  the  hjdrolo^lc 
cycle  but  does  not  attempt  to  solve  a^y  related  practical  problsriSe 

Dynamic  hydrology  is  the  term  for  hj'drologic  studies  involvirg  the 
dynamic  wave  theories  of  fluid  flow ,    These  wave  theories  have  been 
applied  to  overlarid  flow  as  wall  as  oper  chamAel  flow  to  derive  syn- 
thetic  runoff  hydrographs. 

Paraui-otric  hydrology  attempts  to  discover  the  relationships  amo?ig 
physical  parameters  that  are  involved  in  the  partieular  hydrolcgic  even^;s 
that  are  of  interest  az-jd  then  to  use  them  to  sin^alste  non- recorded  events, 
Amorocho  and  Hart  {l$6k )  list  methods  of  correlation  analysis^^  partial 
sysbeiri  synthesis  with  linear  or  nonlinear  analysis;  and  general  system 
synthesis  as  mathods  used  in  parametric  hydrology. 

In  the  prediction  of  peak  runoff  from  a  sf3.11  watershed,  by  a  mathod 
such  as  that  developed  by  the  SCS,  both  stochastic  ana  deterministic 
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hydrology  come  into  play.    The  streamflow  itself,  being  the  result 
of  precipitation,  is  a  stochastic  process.    Watershed  characteristics, 
vhich  modify  or  affect  the  streamflow,  are,  by  themselves,  physical 
paramaters.    Dynamic  hydrology  must  be  considered  in  the  flood  routing 
process  while  parametric  hydrology  is  utilized  to  find  correlations 
among  the  various  watershed  parameters, 

A  discussion  of  the  technique  for  predicting  flood  frequencies 
which  is  currently  in  use  by  the  Soil  Conservation  Service,  is  reserved 
for  Chapter  III,  because  it  is  the  basis  for  this  report. 


■  ,  Chapter  III 

SOIL  CODJSEHVATION  SERVICE  RUNOFF  PREDICTION  METHOD 
Tii3  Soil  Conservation  Sarvica  method  of  predicting  ruroff  from  un- 
gageri  ^jatorsbcds  is  characterised  by  the  devGlopiiiQnt  of  a  syiithetic  unit 
hydrograplie    Undar  this  ciathod  the  Btorm  charactcri^^tics  are  transf criiiGci 
into  a  Bynth'ssised  flood  runoff  hydro-^raph  by  tha  baoin  characteristics^ 
the  shape  of  the  synthetic  unit  hydro^raph  and  the  baseflou  hydrography 
The  general  procedure  \?hich- is  f ollowsd  in  developing  a  flood  rimoff 
hydrograph  is  sho>?ii  hy  the.  block  diagra;;i  in  Figure  lo 

As  sho\m  in  Figure  1^  the  storii  rur-off  vol.r.sa  mvat  first  be  derived 
from  the  stora  characteristics  a^d  tha  basin  characteristics  *    Tnis  storm 
runoff  volirsa  is  then  shaped  irto  a  storin  hydrograph  by  use  of  a  syn- 
thetic unit  hydrograph  laiich  in  turri  is  derived  from  the  basin  and  storm 
characterist ics 6    A  baseflov?  hydrograrh  then  is  added  to  the  storm  hy- 
drograph  as  indicated  in  Figure  1  to  produce  the  flood  runoff  hydrography 
This  sequence  of  operations  involved  in  producing  the  flood  runoff  by- 
drograph  conprises  a  liathsinatical  watershed  inodel.    The  general  frame - 
vork  of  the  watershed  raodel  has  been  incorporated  into  a  program  for 
solu'cion  by  a  digital  computer «     Tha  basin  characteristics  of  a  parti- 
cular watershed  raust  bs  supplied  as  input  to  the  computer  to  transform 
the  general  model  into  a  specific  model  for  the  watershed  in  question. 
Proposed  Uss  of  SCS  ICethod 

The  proposed  way  of  testing  the  SCS  method*^  is  to  use  it  on  an  actual 

■5'«This  ri-3thod  is  outlined  in  the  Ifydrology  section,  of  the  SCS  National 
Engineering  Handbook  ( 196^0 « 
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storm  where  the  input  (stona)  and  the  output  (runoff  hydrograph)  are 
known.    The  procedure  should  be  to  supply  data  characterizing  the  water- 
shed and  storm  characteristics  as  input  to  the  computer  program  which  then 
uses  the  watershed  model  to  simulate  a  runoff  hydrograph.    The  volume  of 
the  simulated  hydrograph  is  compared  with  that  of  the  actual  event.  If 
there  are  major  discrepancies  in  the  volumes,  the  values  characterizing 
the  basin  parameters  may  be  altered,  and  the  computer  program  run  again. 
When  the  volumes  are  adjusted  to  be  essentially  the  same,  the  model  is 
considered  to  be  a  valid  representation  of  the  watershed.    The  validity 
of  the  SCS  method  depends  upon  the  ability  of  the  adjusted  model  to 
accurately  reproduce  the  actual  runoff  hydrograph,  including  the  peak 
discharge. 

After  the  model  has  been  adjusted,  it  is  then  possible  to  route 
design  storms  through  the  model  and  thereby  simulate  design  hydrographs. 
Unit  Hydrograph  Theory  " 

The  unit  hydrograph  (UH)  for  a  given  watershed  is  defined  as  the  dis- 
charge -to-time  relationship  that  yields  one  inch  of  runoff  from  a  storm 
of  given  duration  over  the  entire  watershed  area.    Unit  hydrographs  are 
generally  derived  from  as  many  actual  recorded  storms  as  possible  and 
then  generally  are  assumed  to  represent  the    unit  runoff  functions  for 
all  storms  of  similar  durations,    Sherman  (1932)  is  given  credit  for  for- 
mulation of  the  unit  hydrograph  theory  while  Snyder  (1938)  first  intro- 
duced a  method  for  constructing  a  synthetic  UH  which  may  be  used  for  the 
study  of  ungaged  watersheds.    Several  other  m.ethods  for  defining  synthetic 
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UII  have  been  develop:;:!  since  Snyder  IncJ.udiii^  thoso  by  Cominons  (l9'!-2)^ 
Mitchell  (19^-8);  onu  Gray  (I961). 
Systein  Syp^^si-s  and  ^5£lysis 

WxQTQ  no  actual  hydrographs  are  availabli-,  for  a  watershed;  the  SCS 
itatbod  of  predicting  runoff  can       claBsified  as  a  msthod  of  general 
system  syntbssiSo     If  an  actual  unit  hydrograph  is  available-;,  on  tue 
other  hand;  ths  SCS  rcsthcd  weald  be  one  of  partial  system  syxitheois  vita 
linear  systora  amlysis, 

Systein  synthesis  according  to  Amorocho  and  Eart  (.196-'!-)  I'd  a.  method 
of  describing  ths  operation  of  a  physical  system  with  a  comb irat ion  of 
coj!ip027.e;its  that  exist  in  tha  system  and  uhose  functions  are  kno;?n  and  pre- 
dictable, 

Systeiu  analysis,  on  the  other  hand_,  is  a  method  by  i?h.ich  the  re- 
latioiishlp  betv.'een  the  input  and  output  to  the  systera  is  established 
iTjatherratically  by  rieasurirg  only  the  properties  of  the  irput  (storifi  data) 
and  output  (runoff)  without  regarding  the  nature  of  the  system. 

An  explanation  of  the  linear  and  nonlinear  properties  of  runoff 
prediction  methods  vrlll  follou  later  in  this  chaiDtere 

On  a  gaged  v/atershed  since  an  actual  UII  Is  available,  it  is  an 
arxalyt-ic  (as  opposed  to  synthetic)  function  since  the  internal  charac- 
teristics of  the  systera  are  not  knc.jn  or  synthesised*     Even  though  the 
actual  hydrograph  is  available ^  however;  synthetic  modifications  are 
necessary.     Interception  of  precipitation  by  vegetation  and  baseflov; 
'characteristics,  for  instance^  i-iust  be  synthesized^    Therefore^  even  on 
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a  gaged  -watGrshedj,  t'ae  SCS  method  is  partiality  synthetic^  and  rc.3.y  ha 
clasDified  as  oae  of  partial  system  syiithosis  uith  linear  analysis „ 

Co.  c\n  unpaged  ^jatersh^d.  where  it  is  nccesBRry  to  use  a  Gyrfi;betic 
I3H  the  entire  system  is  synthetic;  thus  the  classification  of  general 
BynthsGis  Kethod  apijlies  to  the  prediction  of  runoff  hydrcgrax)hB  on  un- 
gaged  watersheds e 

The  system  vfnieh  transforVii3  rairifall  to  runoff  can  be  thoughu  of  as 
made  up  of  three  separate  suhsyst 01*153  \-?hen  the  unit  hydrograph  axrprcach 
is  being  considered^    1)  The  first  subsystem^  i.'hich  creates  a  "rainfall 
excess"  function^  ncdifies  the  total  rainfall  input  to  account  for  in- 
filtration, iDterception^  and  depression  storage «    'I'he  rainfall  excess 
thus  computed  is  the  aicount  of  rainfall  that  is  available  for  ruioff . 
2)  The  Decond  subsystem  creates  a  "storia  runoff"  fuxictioo  by  operating  on 
the  rainfall  excess  functioiio    This  subsystem  makes  use  of  the  basin's 
topographic  charaeterictics o    3)  Fina^lly,  the  third  oubGysten^  vjhich  crem- 
ates the  "flood  runoff"  function  uses  baseflo'.-?  information  to  alter  the 
8 term  runoff  function^    The  comblrcd  effect  of  these  three  subsysteKs  is 
assumed  to  duplicate  the  natural  processes  v.'hieh  occui''  on  the  actual 
watershed  being  studied o    Tae  paragraphs  that  follo^r  describe  in  detail 
these,  three  oubsystesis  as  synthesized  by  the  SCS  raethod. 

^££51' •    ''^^^  rainfall  excess  function  (produced  by  the  fir^t 
subsysteiii)  is  given,  according  to  the  SCS  method,  by  equation  (3)0 


(3) 


v^bere  Q  is  tha  amount  of  rainfall  axceso  iu  inches  over  th«  i^al^arshad^ 
P  is  t]i3  storm  raxiifall  in  inches^  £.M  S  is  defined  as  the  maxiiaun 
potential  difference  betv?ccn  P  and  Q  (hence  the  m^iinum  inf iltrr.tion 
cax:7acity)  at  the  tiine  of  the  stom's  begirningo 

Equation  (3)  in  basec^.  on  a  hypothesis.    If  ti.e  equation  can  be 
sho-wn  to  be  true_p  then  the  hypothesis  can  be  aseu^nsd  to  b2  true^  The 
original  hypothesis  can  bs  stat?;d  as  suchj 


where  G  is  the  actual  retention  cluriiig  a  ctora^  S  is  the  potential 
ina>: Imum  retention^  Q  is  the  direct  runoff  (or  the  actual  runoff ) ,  anl  P 
is  the  total  L-torisi  rainfall  (or  the  potential  riaxiKVun  runoff).     So  vjh^^n 
flood  producing  ston'aB  arc  considered^  it  can  be  said  that: 


Equation  (5 )  cannot  be  proven  inathor^ t ically,  honuver^  it  f ollov?s  from 
equation  (5)  which  enipirlcally  hac  been  found  to  bs  valid. 
Since  Q  s  P  -  Q,-  equation  (5)  can  be  rewritten; 
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(6) 


The  SCS  Hydrology  Handbook  begins  the  davalopmsnt  of  equation  (3)  vdth 


the  stateBisxrt  of  equation  (6)0 

As  P  goes  to  irrlTinity^  ths  ratio  (P  -  Q)/S-?-l  a.u6.  q/F-v.-3.  Kaf.p» 
ing  in  mind  that  P,  Q,  and  S  are  total  voliirr.es  for  a  stor^u   (P  -  Q) 
finaIJ.y  fills  tha  entire  coll  profile  to  its  isnxiiauin  value  S  as  de- 
termined for  the  condition  of  the  ijatersbod  at  the  beginniD-_,  of  a  storm, 
Now,  solving  directly  for  runoff  voluna  from  equation  (6) j 

(7) 

A  further  adjustr^isnt  of  P  can       rado  for  the  initial  abstraction 
Ig  which  reduces  P  by  the  amount  of  intercex/blon^  depression  storage, 
and  infiltration  at  the  storm's  beginning  before  runoff  occurs c  Several 
relationships  can  nou  bs  revised  to  allov.'  for  ths  initial  abstraction 
(in  effect,  this  aiiounts  to  redefining  t-ertis  v^hich  nov;  have  slightly 
different  msanirigs  than  the  correspcndirg  terms  in  equaticus  {h)  through 
(7) 

P  ^.  QtG  -ria  (6) 
S  G  -  P  -  Q  (9) 

G  -  (P  -  J.)  "  Q  (10) 
Equation  (7)  can  also  be  restated? 

0  .  J.Ll„3aji^  (11) 
P  »  la-f  S 

Existij^'g  data  from  several  i^atershads  ^fere  studied  and  it  was  found  that 


i 


Ig  can  be  taken  to  be  equal  to  ,2So    By  substituting  this  value  in 
equation  (11);  equation  (3)  is  finally  derived,, 

The-  SGo  has  prepared  a  graphical  eolution  (runoff  as  a  function  of 
precipitation)  of  equation  (3);  i-jhich  is  shoi/n  in  Appendix  A„    A3  is 
eviciczit  in  Appendix  A^  precip5/bation  is  related  to  runoff  through  a 
family  of  curves _j  each  curve  being  drairn  for  a  particular  curve  nmaber 
(CK)o    The  curve  number s  are  a  function  of  S  and  are  given  by  the  re- 
lationships 

...  S-i-lO 

The  parameter  CH  is  usecl  instead  of  S  as  inbcger  values  of  S  do 
not  produce  curves  that  lend  thcniBelves  to  easy  in"uerpo.latione  The 
curve  numbers,  however,  make  the  P  vs o  Q  plot  easior  to  use„ 

The  ruvioff  equation,  equation  (3),  used  by  the  SC'S  is  a  relationship 
betweexi  P  and  Q  involving  only  one  pararLster,  the  parameter  be ing  S»  S 
can  noy  be  empirically  related  to  as  many  characteristics  as  is  thought 
to  be  necessary The  developers  of  this  raathod  found  that  soil  and  cover 
conditions  had  the  most  effect  on  the  value  of  S  (and  therefore  on  the 
ClOc     Tne  curve  number  CH  is  therefore  toriaed  the  soil  cover  coiuplcx 
number. 

For  various  combinations  of  soil  type  and  vegetal  cover,  soil  cover 
complex  numbers  have  been  developed  empirically^    To.  determine  such  a 
number,  snail  V'atersheds  v?ere  found  uhich  had  on.ly  one  type  of  soil  and 
only  one  type  of  cover  condition*    A  number  of  points  v?ere  plotted,  one 
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point  par  stora^  on  a  graph  of  to-Lal  precipitation  in  inches  versi-s  total 
runoff  in  inches ^,    Trie  CIJ  cirrve  that  best-  fitp  this  plot  is  than  taken  to 
be  the  representative  curve  nuiTibsr  for  that  soil  cover  coLipleXo  "  The  re- 
sults of  all  Buch  plots  iD3.cIe  by  th^  SC3  arc  Giij"a:.arized.  in  Appendix  C» 

In  using  the  SC3  Ksthod,  curve  nurabars  (CIl)  are  detcjrrr:in3d  by  re- 
ferring to  a  table  such  as  that  shot-ni  in  Appendix  C,  and  selectirig  the 
CK  that  corresponds  to  the  aijpropriate  ooil  type  and  vegc-tal  cover  con- 
ditiOTie    Tho  so:l  clacslf icat ion  used  clividec  all  soils  into  four  groups 
(A,        C;  and  D)  according  to  thair  permaabillties »    The  portability  of 
group  A  soils  is  highest  and  of  the  D  soils ,  the  lowest o    The  descrip- 
tions of  the  various  soil  types  are  given  in  Appendix  B«    Musgrave  and 
lioltaD;,  as  reported  by  Cho'-?  (196^ h?v3  quantified  infiltration  rates 
for  the  four  soil  groups.    Mini?nuin  infiltration  rates  for  groups  C^, 

and  A  are  given  by  the  ranges  0  to  O.O5  to  Od^^  Od^  to  0«30_, 

and  Oo30  to  Oj\^  inches  per  hour^  recpectivclye 

A  vegetal  cover  condition  is  described  as  some  natural  or  cuJti- 
vated  conditiono    For  instance  the  cover  iray  be  described  variously  as 
fal.lou^  contoured  snail  gi^ain  in  poor  ccndition_j  pasture  in  fair  con- 
dition_;  or  Hoods  in  good  condition. 

The  effect  of  the  antecedent  niDieturc  condition  [Pd^]  is  accounted 
for  by  adjusting  the  soil  cover  ccriplex  nuiifber  previously  calculatedc. 
The  curve  nuiiil3er  vas  originally  developed  uith  average  AM)  being  assuDi^d, 
Kov?ever^  for  conditions  of  very  dry  and  very  v/et  antecedent  soil  UiOistiu'e 
conditions,  type  I  and  type  III,  respactively^,  are  used.     The  typ-3  II 


„  19  - 

condition  Ib  considered  to  be  the  average  condition.    VJhether  the  AMG 
calls  for  typ«  I  or  type  III  tr&atriient  is  detarninod  by  the  guide  linao 
in  Table  lo    ITaa  limits  on  the  dormaiib  sea&on  apply  to  m}fro7~-jn  z^omid. 
aiyl  no  cnov  cover «     ■    •  ■ 

After  tbs  correct  type  of  AM;  is  detorsiinad,;  there  are  standard 
curves  shoyiiig  the  adjustments  to  be  cmde  on  th?  soil  covo-r  cc?;-iplex 
numbers  as  oho-m  in  Figure  2o    Values  for  the  adjusted  CH  c=^n  he  taken 
at  points  between  lines  shown  for  types       11^  and  III  if  it  is  felt  that 
the  AMC  are  suff ieicnt3.y  yell  kno'v?nc 

■    Storsi  Runoffs    The  second  subsystem  (that  which  produces  the  storifi 
runoff  function)  is  assumed  to  be  a  linear  convolution  by  unit  hydro- 
graph  theory^    The  purpose  of  the  convolution  process  is  essent}ial3.y  to 
convert  a  rainfall  excess  function  of  a  certain  volurf.e  into  a  runoff 
function  ha-zirg  idie  same  voJarme^    The  convolution  integral  dcscribirg 
this  runoff  function  c](t)  is  siven  bys  _,.  . 

t ' 

q(t )  ^    r     u(t  -  r )  i  (     )  cl  'i-  (13) 
il  0 

The  kernel  function  u(t  «  'c')  is  the  instantaneous  unit  hydrograph  (tlie 
unit  hydro^raph  nhich  theoretically  results  from  a  rai^ifall  excess  occur- 
ring instantaneously  on  all  parts  of  a  watershed),  i-i?')  is  "t-he  raiin'^all 
excess  input  function,  and  t'   (the  liniting  tine  f or  r  )  is  a  f miction  of 
tQ,  tQ  being  the  dui'a,tion  time- of  the  input  function  (storuO"    For  t-^t^^ 
t'      t  and  when  t  ^-  ±q^.  t'  tQ. 

'Ihe  basic  convolution  process  is  shown  graphically  in  Figure  3«> 
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Figure  3:     Convolution  of  Rainfall  ixooss  i{r)  i^nd  Instan- 
taneous Unit  I-Ivcroar&ph  u(t-t;') 
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Starting  from  tbs  excess  function  i  ( '^f  )  vorGUs      ^  each  inure!"2jrtal 
voluffia  or  raiiifall  ie  liiului plied  by  uha  instantaneous  UH  (XUIl)  and  ic 
suramed  in  tlis  risjht  tiine  seqi"^snce  to  yield  the  storm  runoff  function 
a (t)  vereue  t.    The  instantaneous  unit  hydrograph  Ghoim  is  that  used  by 
the  SCS  i?hen  clavelopiDg  a  Gyirbhstic  unit  h;/drograph  for  an  ungagsd  v?atcr- 
shad.    Thin  oyritliatie  UH  was  developed  for  tho  SCS  by  Victor  MDckus  as 
reported  in  the  SCS  Fational  Engineering  Handbook  {196^1- }«    "J^hG  triangular 
lUH  shown  with  a  dottGcl  line  is  normally  used^  because  of  its  siraplicity 
for  runoff  calculations  ciade  by  hand c    Iioi-?ev€jr,  when  the  computer  is 
utilissd^  the  more  e;cact  cur\rilinaar  lUH  is  equally  easy  to  apply « 

Tha  rainfall  exceeds  function  i( 'J' )  can  be  used  directly  \;ith  the 
lUn  to  produce  the  runoff  function  q(t)e     Hcivevor,  the  way  in  vhich  value 
of  i(T  )  should       calculated  is  not  obvious  from  the  earlier  devalop- 
ment  of  equation  (3)0    Iha  volme  of  rainfall  excGSO  i  (vhich  is  equi- 
valent to  the  vDluiiie  of  runoff  q)  for  the  tima  increment  dr  is  equal  to 
the  value  Q  given  by  equation  (3)  at  tine         dl'   minus  Q  by  equation  (3) 
at  tia'2  X  o    The  voluna  of  precipitation  P  is  the  total  pracipitation  of 
the  storm  up  to  ti^rs    -irt  d  T:'    and   X  ,  respact  ively «    The  inoreTient  of 
(dQ)p  cau.-od  by  tha  stQr:n  of  d length  niiiot  be  given  as  a  difference  of 
Q^G  at  tinos   X\'  CkX  and  Uhls  is  tru2  since  the  runoff  equation  is 

bailed  on.  the  entire  Ptorm  (the  initial  abstraction  nu&t  be  accounted  for 
only  once ) » 

Tha  bace  time  for  the  curvilinear  US  in  Figure  3  is  taken  to  be  5. 00 
times  th3  tiiBs  to  peak.    Inis  value  of  base  tine       aibitrary  but  doas 
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not  incorporate  appreciable  error  in  the  resv.3.ts  as  tris  co.rrespondir;g 
values  of  discharge  are  veiy  smile    The  peak  flov/  rata  for  the  cur- 
vilinear UK  i&  the  saEia  as  that  for  the  triaiigular  vera  ion  ivhoae  base 
time  (deterininsd  mathcmtically )  is  founl  to  13  2c6'[  tiiiics  the  tiriG  to 
peak  (tp)„    The  ■  triangular  hydrograph  hv.sQ  tine  of  2o67  tp  is  used  as 
it  givcB  the  same  value  of  runoff  volura  as  the  curvijlnf?ar  versjon,  the 
volume  heii'Z  represeirlied  by  the  area  under  tao  curve c.    The  peaii  flov; 
rate  q    is  derived  -uoin^  the  triangular  shapad  UH  arid  iB  given  by  the 
ge RG  ra 1  e  qua t  i  on 

vith  tp  s  d/2  -r  i:.  s  :o/:2  ■^  Oc6  to  (15) 

and  t|-,  r  2 067  tp  for  triangular  unit  hydrographs  (16) 

and  t^      5«0  t,^  for  curvilinear  unit  hydrographs  (17) 

v;hsre  K  is  a  cor'stant  whose  valua  depends  o?i  the  irnits  used;  A  is  the 
area^  Q  is  total  voluiia  of  runoff,  t    in  tine  to  paakc  tv,  is  the  base 
ting,  D  is  storm  duration  time,  t^  is  the  time  of  concontratioHp  anl  L 
is  the  basin  lag  tii^:3e    K  is  equal  to  kSh  Vrhen  Q  i«  in  inches,  A  in 
square  mileG  and  tp,  L  and  D  are  in  houirs.    Peak  flow  is  than  in  cubic 
feet  pisr  BGCord^ 

For  the  exatapls  shown  in  Figure  3;  v/herc  the  convolution  integral 
generally  cannot  be  so.lved  exactly,  actual  practice  dictates  that  the 
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storm  input  be  broken  up  into  sraal-l  increments     r  «    Tha  storiit  duration 
D  of  equation  (15 )  is  taen  equal  to  a  r  »    Equation  (1^  )  tO;:;ather  vith 
either  (l6)  or  (IT)  is  sufficient  to  describe  the  runoi"!'  hydi'0£;raph  fxozix 
a  EubwaterKhed  dus  to  an  increnantal  storm*    The  hydrcgraph  caii  be  given 
the  triar!2ular  or  curvilinear  shapa  as  de-sired^ 

Fl£od  5^iS££C'  third  subsyetem  (that  'kjoicb  produces  ths  flood 

runoff  function)  is  the  simple  procedure  of  adding  the  taGefiow'  hydro- 
graph  to  obtain  the  resulting  flood  hydrographo    Tlis  baseflo-.-?  hydrcgiraph 
may  bo  cbtairsd  by  ary  of  ths  standard  textbook  methods  as  lon^  ar  th-^ 
same  procediu'e  is  used  conoistently* 
SiliHBSi  j^lood  Houtin^ 

Mhan  a  vatersh&d  is  divided  into  Gub^/atershads  for  the  purpose  of 
predicting  rimoff  nore  ac-ciu?ately;  it  is  Docessai^y  to  route  hydrcgraphs 
through  tho  natural  channels ^    The  routed  hydrograph  (cutflox-?  hydro- 
graph  frcai  a  cbarmal)  is  then  in  a  proper  foria  to  bs  added  to  another 
routed  hydrograph  or  to  another  Eub'?atersh€d  runoff  hydrography  after  . 
v?hieh  it  Kjay  bs  routed  through  a  second  reach  of  channele 

The  msthod  chosen  by  the  SC3  for  chairiGl  flood  routing  is  tsrrsd 
the  "convex  nrjthod"  after  the  nana  applied  to  ths  ir.atheKatics.l  T)riuctx)3,e 
upon  ^diich  it  is  based o 

The  convex  method  is  based  on  the  principle  that  v?hen  a  flood  ija?:ses 
throug-h  a  reach  of  natural  channel  there  is  a  tire  interval  such  that 

if  \  ^   Oi,  then      ^  Or,   ^.  0^  (3^83  ) 


if  Iji   ^  0^,  then  I-^   ^  0^  ^  0^  (l8b) 

v/here       and  Oj_  are  the  ratee  of  inflo\'  and  outflow-?  respectively  at  time 
t  »  t^;  and  Op  iB  the  rate  of  outflow?  at  tiine  t  xz  t-,-j-it  e  X.^"  Kotie 
that  I2  ( Inf  .low  rate  at      )  does  not  entsr  into  this  r.etbcd.     If   At  is 
chosen  as  specified  then  I^,  O^^  and  Og  forri  a  "convos:  SGt"«    Frcn  con- 
vex set  theory  v.o,  can  then  urito  the  uorkin^  equation 

0^  .  (1  "  0}  o^-j-Cij_  (19) 

whare  C  is  a  parariater  v?itli  the  rar^e  0-^^C-~lo 

Froiii  the  definition  of  a  convex  set,  it  can  be  said  that  yhcnsver 
any  tv/o  poiijts  x  and  y  are  in  a  convex  set   L^then  every  point  on  the 
line  connsctj.ng  thoin  is  also  in         A  formal  definition  froK  set  theory 
given  by  Charnes  and  Cooper  (196I)  states  that i    a  set  of  points  L  is 
convex  if;  yh3never       y-.^L.j,  then  lu  -1-  (1  -  u)y-c.  L  for  all  O^^^u^l. 

The  limits  on  tha  value  of  C  insure  that  aquation  (.I9)  mscts  the 
conditions  of  the   equations  (l8a)  and-l6.b)c    -It  follo-'vs  that: 

It  can  also  be  shov/n  frcn  .Figure  h  that: 

^2  -  Oi        Ii  »  Oi 


At  •    K  ■ 

v/herv?  K  is  a  paraEatcr  with  time  units  not  yet  having  any  physical  in?an» 


s 
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Figui'e  h  I 


Be3.ationshlrj  tstijoen   At  ana  X 


ingc    From  equations  (20)  ai^d  (21)  it.  can       ssan  tnat : 

At-  CK  (22) 

If  use  of  this  mathorl  is  to  be  inade  on  usgaged  xfatershsds,  some 
empirical  relat lQn3hlv)S  noed  to  bs  developed  for  evaluation  of  C  a-icl 
Kc     Thi'6  lias  been  dor^o  hj  the  SC3  from  Kncyledge  gzlne'iJ.  from  stud^^  of 
gaged  v^atorshed^o     K  was  found  that   equ'itioii  (23 )  cxprcssas  tha  value 


C  .  (23  ) 

V  t  1.7 


for  C  very  veil  except  at  extreme  values  'of  V(-Ocp  and  ^10^0 );  v.'liere 
V  is  tbs  average  stream  velocity  in  feet  per  second.     K  \?as  also  founl  to 
ba  the  travel  time  for  the  reach.     If  K  is  in  hours  ar:icl  L  is  the  reach 
length  in  feet,  then:  . 


It  is  obvious  from  the  empirical  rolationshipo  for  C  and  K  that  as 
the  water  suiface  elevation  changes ^  varying  velocities  \vill  cause  C  axid 
K  to  vfi>ry  with  tii^i?)*    A  procedure  for  deterirdnizig  V  '/as  G0uj,ht  go  that 
C  and  K  could  be  assurr.ed  constant  for  any  given  routiw^  prob3.ea.     It  vias 
found  that  if  the  average  velocity  T  vas  based  on  the  stream  flow  velo- 
cities for  all  d:ischarges  of  the  iDfloi-?  hydrograph  greater  thaia  one-half 
of  the  peah  discharge _y  the  resulting  outflov?  hydrograph  vas  nearly  the 
saiiie  as  if  it  were  calculated  using  varying  C  and 
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The  value  of  -At  to  ba  vjsed  is  then  a  product-  of  C  emd  K  and  re- 
Kains  constant  for  the  entire  routings    The  value  of   At  may  turn  out 
to  ba  inconvenient  fractions  of  an  hour  and  therefore  difficult  to  una„ 
SUS  practice  is  to  choose  a  convenient  ^.t '  and  c^.lculate  a  ne\;  C^,  called 
C ,  such  that  routing  by  equation  (19)  gives  identical  resultCc  The 
adjustment  is  . 

C  .  1      (1  -  C)-^-'/^* 

Routing  is  now  accoiT'pliohed  by  retiriting  equation  (19)  to  give  the  resu3,t  i 

0-2  v:  (1  -  C  )  0_j_-;-C'IjL  (26) 

Beservoir  Flood.  Poutinf^  •  .  ..■  " 

•Runoff  hydrographs . are  routed  through  reservoirs  by  the  storage  in- 
dication mathodo    This  mathocl  is  based  on.- the  continuity  equation  shovni 
here  in  its  differential  form. 

i(t)  s  q(t)  (2?) 

v-fhere  i(t)  Is  the  rate  of  the  inflow  at  tiiae  t^  q(t)  io  the  rate  of  out- 
flow at  tiEie  t,  and  dS/jt  is  the  rate  of  change  in  storage*    The  basic 
working  equation  is  yritten;  • 

(lx-i~T2)  ^.tA?  -   {Oy\-  0^)  At/2  ^         "  Sj_  (26) 

vvhere  I  is  inflon,  0  is  outflow; S  is  storage,  t  is  a  time  increment, 
subscripts  1  and  2  are  for  t  -  t-j  and  t      tj  -j-  At  s  tg  respectively. 
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Reduction  of  this  gives  j 
D2finlr:g  connteintB  Cj|  and       eis  i 

%  -    (^i+^'2   "  0;i)At-r2S^  (30) 
.  02Ati.2E2  (31) 
It  thon  follo\;3  that  from  equation  (29)^ 

Q,    .    Cr^  -  (32) 

Cj,  is  first  evaluated  from  known  valutas,  then  set  equal  to 
This  value  of  C^^  is  then  GUustitiited  into  equation  (31 There  is  one 
set  of  values  (O2  and  8p )  on  the  storage -discharge  curve  for  the  reser- 
voir that  saticfies  thir;  linear  eouation»    This  procedui'e  constitutes 
one  £tep  in  solving  the  continuity  equation  by  nuii.^rical  approximt ion. 
To  route  a  hydrcgraph  through  a  reservoir  thi?3  process  muGt  be  repeated 
until  flood  flo..-s  have  nearly  diisjinishod  to  the  original  baseflo'?* 

The  runoff  dynamics  of  an  actual  yatershad  might  be  described  as  a 
distributed^^  non-linear^  time  variant  system,     hany  of  the  models  con- 
structed to  describe  this  system^  horever^  are  classified  as  lu;i'ped, 
linear,  and  tise  invariant ,    Early  models  \iith  these  assuinptions  v?er«3 
those  developed  by  Sherimn  (1932)  and  Bernard  (1935  )o     Sheman  developed 
the  unit»graj)h  concepts  often  callad  unit  bydrcgraph  theory.    Bernard  is 
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given  the  credit  for  ;L5npleEant ing  xxbc  of  the  dietrlbui; ion  graph  (a  bar 
graph  C0::2parable  to  a  imlt  ]:ii'drograph  except  that  runoff  voluiaes  are 
expressed  as  porcom;  of  total  riuioff  versus  tiiia),, 

A  distributed  systeia,  in  the  field  of  hydrology,  is  defined  as  one 
with  time  and  space  variation  in  rainfall  intensities e    A  luiapc;d  systi?i\ 
or  model  if:  one  vhieh  allows  only  time  variation  but  no  space  variation 
in  rainfall  intensity.    This  is  to  say  that  rainfall  intensity  is  uniform 
over  the  v;atershed  area  unasr  eonsidaration  at  any  c:iven  tij;;ee 

The  SCS  msthod  in  its  basic  form  (as  currently  explained)  yould  he 
classified  as  a  luiiipsd  systeuio    By  using  the  SCS  runoff  prediction  mat  hod 
in  conjnnctioa  uith  reservoir  B.n6.  cha-niel  flood  roircing  techniques,  the 
resulting  systen  x-rill  in  part  be  distribntedo    Cho-.'  (196"(')  describes  this 
type  of  system  as  a  "distributed  systcji  of  lunipcd-system  models "c  On 
each  of  several  subv^atershods  the  rainfall  intensity  is  uniforra,  or  each 
of  the  subsystems  is  lumped.     However,  at  any  given  tvi^ia  ^  the  intensity 
can  vary  from  subwatershed  to  subwatershed  and  is  in  this  sense  distrl° 
butedo 

A  system  is  linear  if  the  miit  responses  of  various  in^^ut  rsagnitudes 
are  identical*    The  ^^atershsd  runoff  function  is  linear  if  one  unit  hydro- 
graph  is  sufficient  to  accurately  describe  all  storms  of  a  given  duration 
occurring  on  a  given  t;atershedo    This  is  never  the  case,  ho:?ever,  as 
Chills  (I95S)  and  others  have  sho;?n<,    Childs  states  that  the  larger  the 
flood,  the  higher  the  peak  and  the  shorter  the  tirae  to  psali  on  the  unit 
hydrographo    He  derronst rates  this  fact  by  showing  four  unit  hydrcgraphs 
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of  ubrGe'-hour  siorvn  on  I^,U3atuck  Kivur  at  Thomnton;  Conv^octicirb .> 
The  peak  of  a  iiiajor  flood  has  more  than  double  the  peak  6.^ri.Ye&  for  a 
minor  flood  aad  about  t^o-tliirds  of  the  tirio  to  peak.     Obviouoly,  tha 
variance  of  the  storm  pattern  or  raixifall  excess  function  affects  tha 
unit  hydrograph  shape  o    Kon-linearity  is  then  the  case  vjhen  the  clKrivGd 
unit  hydrcgraph  varies  \jxtb,  ctorid  pattorii  and  tcrbai  rainfall  excess  araounto 

A  systsm  that  gives  the  Gam3  reFpDiiso  for  a  giYGn  inp'it  recardlosB 
of  tiBio  is  a  time  iDvariaiit  system,    Tho  Gub^-?atershed  portion  of  the  SCS 
modol  is  based  o;i  unit  hydroaraph  thaory  and  is  th&refore  ti&e  invariant . 
This  is  not  too  restrictive  an  assuKption  as  the  rair-fall  inpixb  to  tlie 
subyatersbad  runoff  mode 3.  is  preadjusted  for  soil  ani  cover  conditions, 
antecedent  moisture,  and  time  of  year^    Those  arc  the  niain  iteiiis  giving 
the  ruDOff  function  {system  responye  )  a  tiria  variaiit  nature. 

Spatial  and  teEporal  variations  in  rainfall  disti ib^rtion  can  be 
acccuiited  for  vn.tb  the  SC3  approacho    Thsse  variations  are  simulated  by 
dividing  a  v/atershed  into  SGveral  cubw'atersheds  and  by  r^akin^  the  time 
incrvinont  on  tbe  precipitation  versuG  tiEio  input  graph  snail  enough  co 
that  the  graph  nearly  Eatchos  the  actual  curved  relat ionshipo    The  total 
amount  of  rainfall  can  be  different  for  each  subyatershad.    The  storm 
pattern  can  have  an  infinite  number  of  shape  variat ioiis .    Spatial  var- 
iation 18  not  total3,y  accounted  for  since  there  is  still  no  spatial  var- 
iation allovjed  y5.tbin  the  subyaterfihed;  thus  cpatial  variation  is  only 
partially  ciriulated.    The  subyatersheds  can,  boy  ever,  be  Fade  as  small 
as  is  practical  for  the  time  and  expanse  that  can  be  spent  in  the  analysis « 


The  tine  of  storr-  besiBniD3  and  fiu::'at5.on  can  also  be  cleslgr'ated  ac 
necessary  to  conrorm  to  an  actml  ctormo  Thsrafore,  teiapcral  variation 
is  also  partially  slrjiAlateda 


Chapter  lY 

SOIL  CONSETiVATIOH  SERVICE  PROJECT  FOmmjh'IKM  PPvCCiBAM  »  BYDHOIOGY 
The  Soil  Conservation  Service ,  in  1963;  obtainerl  (from  a  private 
consulting  f  inu  Knov/n  as  C«  EU       R.,  Inc.  )  a  digital  computer  pro^jrain 
deBignad  to  solve  the  "runoff  hydrograph  preclictiGn"  relationships  df;s- 
cribed  in  Chapter  III.    Tha  computer  prograa^  uliich  \vas  writtsn  in  the 
Fortran  language  for  the  IBM  7090  computer  system,  is  entitled  ''Project 
Formulation  Program  <»  Bydrology"  and  ie  a  i7i3.thnmtical  model  of  a  ^-^ater 
Bhed  ejiparienclng  p2:'«cipitation,  Burface  runoff,  and  Cuann^l  flou^  Ths 
program  can  rouce  an  unlisiitcd  number  of  hypothetical  or  actual  storns 
through  a  vatarshed  having  as  mny  as  6o  structures  and  120  cross-sec- 
tions described. 

The  program  has  gona  through  several  revisions,  most  recently  in 
October,  1957 .    It  ic  now  written  in  J'ortran  IV- G  and  utilizes  the  IBM 
360/40  coiTiputer  system.    The  most  recent  instructions  for  using  the 
Project  Fcrraulation  Program  (PFP)  uere  published  by  C,  E,  I.  R»,  Inc, 
(I96H)  and  the  SOS  (I965)e 

The  PFP  consists  essentially  of  a  series  of  subroutines  v/hich  Esiy 
be  solved  in  a  variety  of  combinations  and  sequences;  a  schenatic  over- 
view of  the  PFP  is  slioun  in  Figure  5»    The  pa.rticular  sequence  to  be 
followed  in  a  given  problem  is  detarminad  by  instructions  to  a  "StaU" 
dard  Control"  routine,  uhich  calls  the  subroutines  in  the  proper  order, 
and  provides  the  additional  inforiTation  necessary  to  pc-rforr-i  the  york 
of  the  subrouline,    0?he  entire  program  is  monitored  by  an  "Executive 
Control",  vjhich  describes  the  particular  storis  under  consid?ration. 
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To  foriimlate  a  PFP  problo!!i,  it  is  first  moss^ary  to  read  in 
several  ijable&,  as  indicated  on  Figure  5^  soria  general  and  sous  spe- 
cific for  tha  par'fcicuj-^r  i-?ater6hsd«    After  this  ic  done,  Standard  Con- 
trol cards  arc  read  in  to  describs  the  logical  cequence  in  tmich  flood 
routings  thrciigb  the  stream  reaches,  and  ctructuxes  should  bs  perfornodo 
This  sequorce  of  opsrations  is,  of  coiu'-se,  idontical  .to  tiio  sequence 
that  the  hydrologlst  follov?s  in  solving  the  problem  mani^ally^ 

To  reflect  a  future  charge  in  a  i/atershed,  any  portion  of  tha 
Standard  Control  can  be  modified  at  any  tine  by  addition,  altoration,  or 
delation  of  v'atershed  data^    Executive  Control  cards  are  then  road  imo 
tha  computer  to  des cribs  each  of  the  altcrrative  storia  situations  that 
are  to  be  analyzed »    An  infinite  nnsiber  of  problems  niay  ba  sat  up  and  run 
by  the  Executive  Control  for  a  given  set  of  ^^atarshad  characteristics. 
Each  directive  of  the  Execul^ive  Coxitrol  spacifies  tha  stori'i  precipi- 
tation, its  starting  tima,  the  antecedent  Koisturc  condition  under 
which  it  is  to  be  analysed,  and  the  portion  of  ths  watarshad  through  which 
it  is  to  be  routed = 

Each  subnatershed  described  may  have  different  total  rainfall  amounts, 
different  storm  start i?]g  tii^ss,  different  storm  tima  durations  and  pat- 
terns,  and  varying  antecedent  moisture  conditions,    A  given  storm  sit- 
uation is  than  applied  to  tha  Standard  Control  to  effect  a  solutiono 
Additional  storr^s  ray  then  be  set  up  and  solved*    As  has  been  seen,  tlia 
Standard  Control  can  at  any  time  be  changed  to  reflect  proposed  changes 
on  the  watershsd.    After  a  change,  any  hypothetical  or  actual  storm  can 
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again  bs  applied  to  the  modified  imtcrshGcl  to  evaluate  ths  changaG  nacia 
in  the  oiitflQi^  bvdrcgriipli. 

The  SCS  program  has  a  total  of  twelve  subroutines  (seven  of  which 
are  computa t ional ) ,    Their  relatioiishlps  to  the  main  program  arl  to 
other  subroutines  are  shoifn  in  Figure  60 

The  seven  compiit  at  ional  sybroutlnes  are  called  RilOFP'Z.  EKSVQR^ 
REACH,  hmim,  OUTPUT,  XNTEBP,  and  PEAK.    They,  respactxvely,  calculato 
a  -runoff  hydro'^raph  for  a  subv;atershed|  route  an  input  hyai'ograph  through 
a  reservoir  by  the  stora5,e-indicatiOii  mathol;  route  an  inpiii;  hydrograph 
through  a  channel  roach  by  the  convex  method;  add  tvo  hydrographs ;  print 
and/or  punch  tba  output  options  specified;  adjust  the  ti^ic  increraent  for 
RESVOPv,  REACH,  ard  ADIWD;  and  calculate  coordinates  of  the  highest  peaks 
(to  a  raxiiTiUin  of  ten)  of  a  hydrog.raph. 

The  other  five  subroutines  are  called  SAYKOV,  UPfATE,  SETUP,  READ  5, 
and  READIUo    SAVKOV  gives  the  writer  of  Standard  Cent r el  the  ability  to 
move  an  entire  hydrcgraph  from  one  .storage  location  to  another  and  to 
save  it  there  until  needed  later*    UPDATE  causes  the  Standard  Control 
lietiug  currently  In  force  to  be  printed  out,    TniB  woulri  ba  used  after 
Eome  alteration  of  Standard  Control  had  been  loadee    SETUP  ha>:;  the  over- 
all control  of  the  coifipu^jat ional  nubroutineSc    READ  5  ^^^^^  READIN  cause 
the  program  to  read  in  the  data  which  describe  the  physical  features  of 
the  ■.?aterehed  and  the  properties  of  the  storm  being  simulated. 
Tabular  ^J-ta 

Several  types  of  data  that  are  necessary  for  the  operation  of  the 
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Project  Formi'lation  Program  can       suppilod  to  th::  eompiruer  in  tabular 
forme    These  data  ara  vcluas  for  tha  relationship  b£;tv?eon  the  rcutirg 
ccefriclQiit  C  and  avarage  velocity  Vj  the  curvilinear  unit  hydrographj 
stcrm  patterns  J  the  gt age -g t  orage -die  charge  relation  for  reservoirs  and. 
tba  stage-area -cliGcbarge'  relation  for  charjiel  cross-scctionG , 

The  C  -  V  graph  is  tabulated  for  - corr-puter  use  in  velocity  injre- 
niants  of  0»2  feet  per  SGConls    Thi^^  unit  hydrograpb  is  supplied  in  dimen- 
sionleaa  form  a3  q/q^^  vorsuG  t/t.^c     Design  storm  patterns  are  der^cribsd 
in  cumulative^  d line ns ionle s s  form  with  coordinates  of  P/Prp  and  t/tp 
vhere    psp  is  the  storm's  total  precipitat lone    Actual  storms  can  be 
read  in  either  dimensioned  or  dinians ionle S:2  fcrni  but  alv/ays  an  a  cumu- 
lative grapho     If  actual  dimensions  are  used^  rainfall  depth  and  dnraticn 
in  Executive  Control  are  to  be  given  as  loOo     If^  ho^'^ever^  dimensionless 
coordinates  are  used  for  actual  storms  in  tabular  data,  the  actual  values 
of  do-pth  and  duration  irj.st  be  supplied  in  the  Executive  Control,  Tne 
structure  and  cross-section  data  can  be  given  in  any  desirable  increii^airb 
of  elevation.     Tills  should  be  a  SKall  enough  increiiiant  so  that  the  assumed 
straight  line  between  points  is  nearly  coincident  v?ith  the  curved  re- 
lationship from  nhich  the  points  were  selected.     The  elevation  incre- 
ment  can  be  varied  et  i/ill  throughout  these  tables « 

The  vay  in  \?hich  these  tabular  data  are  utilized  by  the  Executive 
and  Standard  Control  is  outlined  by  the  flow  chart  in  Figure  5c  The 
probleps  posed  by  the  Executive  Control  arc  individual,  distinct  storias, 
•an  infinite  nuxnber  of  v?hich  can  be  processed  by  the  ??P*    After  problem 
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KOo  1  has  been  vjorked  by  the  Standard  Control  tbroush  Its  various  sub- 
routinas^  problam  No*  2  is  tbc;o  described  by  the  Executive  Control  airi 
solvad^    This  contiviueR  until  all  stated  problems  h^.ve  teen  BolveO.^ 
These  subroutines  in  turn  use  the  tabular  data  (shoun  in  squares  o'l 
Figure  5)  to  conplete  their  function,    REACH  needs  the  C  -  V  and  the 
cross -sect ion  graphs  to  ccirplete  its  task .    The  croES-:jecti03.i  iuf orrration 
is  needed  uhether  OUTPuj?  is  called  or  rote     RIICFFX  Uf?es  the  UII  and  either 
a  design  or  actual  Rtori.n  P  -  t  graph*    SAVIIOV  and  APDillD  do  not  rake  v.3e 
of  the  tabular  data„    KS3VCVR  must  call  the  structure  data  to  perform  ite 
routing  problems    All  five  of  the  working  subroutines  siay  have  to  call 
OUTPUT  (if  specified  by  standard  control)  before  returnii)^;  to  Standard 
Control c 

Standard  Conti'ol  is  an  ordered  sot  of  instructions  to  the  coiaputer 
that  outline  the  sequence  of  operations  that  simulate  runoff  from  a  \?ater" 
shed  areac    These  operations,  which  have  already  been  outlined,  are  done 
by  the  subroutines  RHOFFX,  RSSVOR,  REACH,  ADDHYD,  and  SAVMOV.    Figure  7 
sho'ws  a  sample  watershed,  with  two  reservoirs  and  two  main  stream  branches 
which  is  subdivided  into  five  subwatersheds .    A  table  of  the  correspond- 
ing Standard  Control  listing  for  this  sample  watershed  is  slso  givsn^ 
The  structure  and  cross-section  locations  are  numbered  so  that  the  reader 
can  follOw  through  the  Standard  Control  to  see  that  it  adequately  describe 
the  way  in  which  rainfall  excess  is  converted  into  runoff.    There  are, 
however,  sone  assuinptions  in  this  connection  i/hich  make  the  simulation 
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Figure  "J  i    Baijplc:  V/atershad  with  CorrespoiidiDg  Standard  Control  for  the 
Soil  Ccncei-vat ion  Service  Projtict  Formulation  Prograiu 
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much  easier  v;lthoirG  acid  ins  subotaiit  ial  eiTor  to  the  ops  ration.    For  ' 
example,  tha  outflou  hyaro^raph  from  structirre  1  is  routed  tbrough 
area  II  as  if  no  additional  runoff  yere  beiD^  supplied  from  area  IT„ 
Now,  tbare  exists  an  oubfiow  bydi^ograph  from  the  channel  reach  in  area 
II  at  section  1,    To  this  b^drograph  is  then  added  the  runoff  hydrograpa 
from  area  II  ao  if  it  bad  not  arrived  at  section  1  by  vay  of  the  channel 
BysteE\»    Tliis  assur;:j:>t ion  ie^  of  course^  not  phy;3tcally  the  case  but  soon 
after  leaving  ths  upstream  subnatersheds  the  channel  hydrograph  is  very 
large  coinparcd  to  the  local  area  runoff  so  that  the  error  retiulting  froDi 
this  asEU;;ri)tion  is  considdred  snails, 

Computer  locations  for  bydro^raph  storage  are  eaonn  on  the  table 
included  in  Figure  7«    There  arc  seven  different  locations  supplied  \-?here 
all  the  coordiratejj  of  points  from  a  hydrograph  can  bs  stored  (up  to  500 
points ) .    Those  points  are  all  described  by  a  1 5-rr.3  and  discharge  value . 
All  subroutines  have  an  inpub  and  an  output  hydrograph  except  ADDIiYD  \vhich 
has  t'.?o  input  and  ona  output  hydrographs ,    The  nunibcring  of  storage  lo- 
cations as  shOi-rn  is  set  and  consistently  used  except  for  the  SAVMOV  sub- 
routine vhiich  can  then  be  used  to  place  a  hydrograph  in  the  correct  3.oca- 
tion  for  the  tiibsequenb  operation.    This  rigid  nuuuisring  system  is  used 
GO  that  fcuer  ei'rors  will  be  iiiade  by  a  technician  pieparing  the  Standard 
Controls 

A13.  the  possible  output  options  are  also  shonn  in  Figure  7»  Thsso 
include  printing  out  all  the  values  of  the  hydrography  the  corresponding 
v?ater  surface  elevations,  up  to  the  ten  highest  psah  dischargs  values  and 


their  times  of  occurrence^  and . tha  total  voluir.Q  under  the  hydrogreph 
givrif  In  tb:>;'eG  oeto  of  urs^tD;;  cfc-hrG^  acro-feet^  arci  inches  of  r\ui,off 
over  tha  vjater^^heae 

Cthsr  irforEatioH  then  that  Ghov;n  must  alno  ba  supplied.  Area 
in  square  railes^  runoff  curve  number  CK_j  and  tins  of  cones ntrat ion  in 
hours  imist  also  be  given  yith  tha  ftHOFFX  coiicnando    RESYOH  Dast  also  kno''; 
the  surface  elevation  of  ths  re&ervoir  at  tim-^  t  »  0.     lc-r)3,th  of  a  reach 
in  fcst  and  a  routing  coefficiant  (if  ona  1g  to  be  supplied  insteuci  of 
being  calculatad  by  tn;  progrom)  are  n';eded  to  eoriplete  the  EEACK  sub- 
routinoe    ADDHYD  and  SAVMOy  are  coicplets  as  shonn  in  Figure  7^ 

With  regard  to  the  numljering  of  stream  channel  cross "Sacticns  anc*. 
control  Gtructiu'es,  numsrical  valuea  in?.y  be  given  in  any  ordor^  althou;3h 
a  logical  sequential  ord^r  is  easier  for  others  to  follcn,    Ho^/ever^  the 
uay  in  which  these  mmiborb  are  used  in  Standard  Control  is  very  itnporbant 
as  will  bo  seen  when  the  use  of  the  Executive  Conorol  is  Gxp3.ainedo  Tlie 
iiiportant  thing  to  rei:eniber  is  that  any  number  used  more  than  once  must 
appear  consecutively  in  the  Stan'ilard  Control  listiiigc    That  is^  if  crcf!3- 
section  001  is  used  a a  it  is  in  Figure  7,  the  four  tines  that  it  appears 
imst  follou  one  anothar  without  interruption  by  other  nu;:Vbers<, 
Executive  Control 

Executive  Control  is  a  oet  of  coinraands  that  have  the  over-all  con- 
trol of  the  computer  solution  to  stated  hydrologt'C  problems .    As  in- 
dicated earlier _p  each  E;cecutive  Control  directive  specifies  the  proper- 
ties of  the  hypothetical  or  actual  storm  that  is  to  be  rontedc  Stor;^i 
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properties  accounted  for  in  this  model  incluas  storm  pattGrn-^  total  ralji™ 
fall  amount^  stem  length,  and  startiug  tim^o    IfnQ  propsrtic?G  x:.d.y  varj 
from  subv/atershed  to  cv!b":atershecl  but  are  assunsd  uniforn  ever  eacli  in- 
dividual subwatersbed. 

The  storm  pattern  is  genei^ally  read  into  the  computer  v?ith  a  dircen- 
6ionlest>  curaulative  rainfall  amount  ordinate  and  a  dimajisioned  tims  co- 
■ordinate.    However,  it  is  possible  to  read  in  both  coordinates  yith  either 
dimenBia:;9d  or  d ime ns ion].es s  scales,    ^he  dinsriDionless  valui^s  are  then 
traiisformsd  to  dineusioned  values  by  specifying  the  total  rainfall  amoum}. 
the  total  storm  leir;th  tbat  applies  to  the  particular  case,  or  both.  This 
flexible  arrarigemaDt  is  useful  for  applying  hypothetical  storr:3  to  the 
system^     In  the  case  of  actual  storrs,  cmuulativc  graphs  vrith  the  correct 
time  and  precipitation  dimensions  can  be  read  in  directly, 

II^REM,  COi-D?lir,  ErlDCl'IP  are  additional  statements  used  in  Executive 
Control  which  simply  are  shorthand  for  incrcmrirt,  compute,  and  end  of 
computations,  respectively.    They  are  not  subroutines  but  simply  commands 
that  are  used  by  the  main  X)rogram  to  carry  out  the  desired  control  se- 
quences, 

Examp^les  of  |dL££}'^.^i!l£  £2ii'li£2P"*    Several  example  problemis  are.  il- 
lustrated in  Table  2  to  demonstrate  the  flexibility  of  the  Executive  Con- 
trol,   Problem  Ko,  1  breaks  the  watershed  into  three  subuatersheds  each 
receiving  different  amounts  of  rainfallo    hainfall  Table  ho.  1  Hould  be 
given  in  dimensionless  depth  units  and  time  in  hours  since  th:;  amount  of 
'rainfall  is  specified  and  1.00  is  given  under  rainfall  duration.  The 
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antGcodent  Koi store  coiiditlon  (All!)  is  type  II  or  nonral.    Starting  time 
is  given  as  zero  on  all  subuaterebeds 3 

Problem  No,  2  is  Bimilar  to  NOo  1  except  that  storm  movement  is  - 
accoui;ri;eci  for  \yj  varjins  the  start  ins  tina  on  the  various  suu^jatershects » 
AlGC  tha  storm  pattern  is  considerod  to  ba  different  in  shape  cr  lergtb 
ao  rainfall  Table  No,  3  is  used.    The  AMC  is  typo  I  or  very  dry. 

Problem  No«  3  is  an  actual  storm  that  is  taken  as  unifoi-m  in  -ainount 
over  the  entire  vjaterslisd,    Tcs  starting  tiroa  varies  for  two  ii^ubareas  ana 
a  new  lable  (Ho,  2 )  describes  the  s terra  patterns  ^vhien  is  given  in  true 
units  of  incbes  and  hours  since  the  factors  used  for  rainfall  depth  and 
duration  are  both  given  as  unity^,     Pi'oblens  k  and  5  ^2;'e  design  storms  and 
are  uniforra  ever  the  entire  area  in  rainfall  aiaount  and  stare ing  t ins » 
The  rainfall  depths  given  would  correspond  to  so2\€  calculated  design  s terra 
amount  for  a  given  frequsncy  storm  for  tbe  area  in  quest iono  Rainfall 
TablGS  k  and  5  niig^it  be  storri  patterns  for  12  hours  and  2h  hour  storr.iSe 
A  nomal  A  MO  would  aluays  be  assumsd  for  design  storr-is, 

IHCREM  states  the  time  increm^nb  in  hours  that  the  coseuter  should 
use  in  developing  the  hydrotsraphs ,  COIvPUl  instructs  the  cc:nputsr  to 
coEipute  what  is  called  for  in  Standard  Conbrol  bstyecn  the  structuro(s) 
and/or  croGS-sectlon(s  )  indicated  in  the  command.  EliDCl'lP  indicates  the 
end  of  a  COMPUT  series  and  sends  control  back  to  the  bagiiiaing  of  Standard 
Control  to  starb  a  nev-  runoff  prcibleii;  based  on  anotber  storm  described 
in  the  Executive  Control,," 


•    Chapter  V 
FRSPARATION  OF  WATfiRSKED  DATA 

ThQ  preparation  of  vjatershed  and  storm  data  for  us  a  uith  the  Soil 
Conservation  Service  Project  Formulation  Program  (PFP)  will  nov?  be 
deBion^^trated^     The  example  used  will  bs  the  Diick  Creek  IfetGrobed  near 
Erockway_j  Montana  and  the  ntorm  of  June  3^6 ^  jSS'}^    Fic;ura  8  is  a  ran 
of  the  Duck"  Creek  dralrage  shoi/ing  the  location  of  rain  gages  and  v/ater 
stage  recoi^derSc    The  streaya  nctuork  is  also  clearly  seen^,  ac  is  ths  sub- 
division,  of  the  v?ater3lu^d  into  sub'./atershsds , 

The  largest  of  the  four  v?atercheds  under  study  by  the  Drainage 
Correlation  Research  'Project^  I>jck  Creek  (54  cquare  miles),  located  south 
of  Brockv;ay,  Montana^  v;as  selected  for  the  investigation  reported  liereino 
It  wa8  chosen  for  several  reaBons  ;    1)  its  £ir:'.e;  2)  the  occurrence  in 
1965  of  a  rain-caused  runoff  event  which  appeared  to  be  Buitable  for 
analysis;  and  3)  particular  interest  in  tliis  uaterchcci  by  the  Montana 
State  office  of  the  SC3.     Instruments  installed  at  Duck  Creek  by  the 
Drainage  Correlation  Research  Project  include  three  i/eather  stations 
which  rnaasure  mIvjI  velocity^  air  and  noil  temperature ana  soil  noistiir-^ej 
four  non-recordiug  and  three  continuon3-r6Cordir.g  rain  gageoj  and  two 
continuous -recording  vatcr  level  recorders.     One  uater  level  recorder 
located  on  the  East  Fork  of  Duck  Creek  uhsre  it  is  crossed  by  State  High- 
uay  253.    The  other  is  located  on  Duck  Creek  530  feet  n03?th  of  the  center- 
line  of  the  county  road  between  cections.  11  and  ih ^  T17H,  RU6£c 

The  PFP  ie  a  generalised  program  designed  to  be  used  on  a  wide 
variety  of  vmterebeds,  and  ie  capable  of  processing  many  types  of  data . 


Analysis  of  the  Diick  Creek  watershed  did  not  require  utilization  of  £-.ll 
the  program  capiibilit ieo ^     The  data  supplied  to  the  computer  for  Duck 
Creek  relate  to  1 )  cubwatershed  areas  and  cl'iarnel  leDgths^  2)  etago- 
area  and  stage -discharge  ralat ionshl ps ,  3)  curve  niuubers  and  antecedent 
nioisture  conditions,  'i- )  ti;iie  of  concentration,  and  5)  storm,  data,,  These 
are  dlfrcussed  in  detail  h'alo\u   :     h    '  '  •     ■  '■  ■  -     ■■     •     ■   -  ^ 

§H^iif^^5£^J±£^  -5±3f-'£         Channel  La iK^-hs 

For  purpoGS5  of  constructing  a  model  of  the  Duck  Creak  channel  net- 
work;  the  Duck  Creak  uatershed  uas  subdivided  into  tv?cnty»f our  sub-fater- 
shedo  o  . 

The  projected  area,  channel  meander  length,  and  the  flood  plain 
length  Mere  neasured  for  each  of  the  subi^atersheds «    Arear  were  determined 
by  planiiTieteriiig  aerial  photographs  (scale:    ^-i-"      1  mi=)  supplied  by  the 
SOS.    Both  meander  and  flood  plain  lengths  were  measured  on  the  same 
photographs  with  a  map  and  plan  measure^ 

For  the  outflow  section  of  each  of  the  subwater-sheds ,  stacre-area  and 
stage "discharge  curves  muot  be  supplied  to  the  PFP,    To  get  this  infor- 
mation it  was  necessary  to  obtain  channel  cross-sections  in  the  field  at 
twenty-six  locations.  ,  ' 

From  the  eross-section  survey  notes,  stage-area  curves  were  obtained 
by  assuming  that  the  water  surface  is  always  lateralJ,y  horizontals 

The  stage-discharge  curve  for  each  section  was  calculated  by  the 
M^nnirg  equation; 


n 

\vhere  n  is  the  L^aiining  rougbnoss  coefficient,  R  ths  hyuraulic  radius,  and 
S  the  friction  j?lopa  at  that  Beet  ion, 

Fr3.ctiqn  Slops  s    Friction  slopa  is  taken  to  ha  equal  to  tho  cbamjel 
filopa  on  Duck  Creek  for  several  reasons* 

Firrt  of  all,  Eonderson  (I966)  etatas  that  because  of  the  steep 
slopes  u&ually  presera,;  in  the  upper  eatelimant  regi0i]:3  yhera  the  runoff 
problem  exists,  bed  slope  v/ill  be  tho  only  significant  slopa  term  ancl 
the  di;i;charge  will  be  a  function  of  depth  alone «    Bsd  slopG  can  ba  con- 
sidered, the  significant  slops  factor  for  slopes  greater  than  10  feet  per 
miJ.e^ 

Fnrtherraore ,  the  cross-Gcct ions  on  D^jme  Creek  are  taken  at  an  aYcrago 
of  ti;o  miles  apart o    Also,  thare  are  many  variations  in  cross-section 
ehape  betv^een  any  two  maasured  cross-cections  c    Therefore,  the  actual 
friction  slope  could  be  expected  to  deviate  from  the  average  bed  sloiJe 
throughout  tlie  channel  reach »    The  diD^'ference  betv7ec?i  the  average  friction 
f^lope  and  the  average  bed  slope,  'nouever,  vjould  be  expected  to  be  very 
Eu'-all.^ 

Zcme  of  Tran.riltJ.cnis    Ihe  stage-discharge  curves  used  are  a  cosibina- 
tion  of  curves  based  on  different  assuiipt ions «    The  channel  slope  \med.  for 
loner  flows,  vhere  the  flow  is  confined  to  the  main  channel,  ic  the  elope 
of  the  channel  bed^    For  very  lai-ge .  f -Ions  covering  the  flood  plain  and 
floning  at  depths  greater  than  the  vegetation  height,  the  effective  s?lope 


is  tbaij  of  the  flooil  plainc    The  inportance  cf  this  con'jiclerat ion  of 
slopa  is  apparent  vYien  the  ratio  of  K'jander  length  to  flood  plain  length 
is  high;^  ap  it  is  on  louer  portions  of  Duck  Creek  vjhera  it  reaches  a 
value  of  t'.:'D, 

A  transition  zone  exists  in  i-?hich  the  effective  slops  varies  ba- 
twe^n  tlie  channsl  bed  and  flood •  plain  slopes'  (shonn  in  Fi^^ure  9)«    This  ■ 
transition  would  ba  difficiU.t  to  sir.iu3.ats  v?ithout  extensive  knov-fled^ce  "of 
the  uatershed.    Therefore    in  this  study  the  Btago-discbar^je  curve  is  a 
fuucticn  of  tha  channel  bi-^d  slops  until  approximately  on5  foot  of  v?at(?;r 
is  flO';^ing  on  the  f j.ood  plain;  above  this  depth  tha  flood  plain  slope  iB 
co.a3 idered  ef f e ct  ive  ^ 

£5H^^£2JH£2  Coefficient  t    The  vaj.ues  of  I'^inning's  roughnesB 
cosfficiGnt  "n'"  vere  ss3.ected  nstrii^  s-  nmiber  of  photographs  of  channel 
sections  for  which  "n"  had  been  calculated  by  Vo  0?.  Chow  (19^0),  These 
were  compared  VJith  photo^jranhs  of  Du.ch  Creek  to  deternrine  the  "n"  value 
used  for  the  various  channel  reaches  of  Duck  Creek, 

An  att:n;pt  was  laade  to  check  this  dsterminsition  with  a  set  of  U= 
Geo.lO{r,ical  Survey  (U3G3)  stereoscopic  color  slides  that  have  been  used 
for  roughness  coefficient  deterriinat ions  for  that  organisation.  IlhTor- 
tunately,  tha  set  of  slidas  viev-ed  did  not  include  exa:riplcs  sirailar  tc 
those  found  on  Duck  Creek. 

One  further  comparison  was  iiiade  in  an  atteiapt  to  evaluate  the 
effective  roughness  of  Duck  Creek,    This  ir-etbod  was  developed  by  Cowan 
(3.956)  and  can  be  suTiiuarised  with  the  following  relationship: 
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n  =  (dq  i-        +        -I-         -f       )n^  (3^1) 

where  n  is  the  desired  roughness  coefficient  to  be  used  in  Manning's 
formula,     uq,  n-^,  112,  n-^,  and  n^  are  the  fractional  values  of  n  deter- 
mined respectively  by  type  of  channel  bed  and  side  slope  material,  effect 
of  surface  irregularities,  variation  in  shape  and  size  of  cross-section, 
effect  of  obstruction,  and  vegetation  conditions,    n^  is  a  multiplying 
correction  factor  and  is  a  function  of  the  ratio  of  meander  length  to  . 
flood  plain  length  for  any  given  reach  of  channel.    The  roughness  coef- 
ficient given  by  this  method  is  larger  than  those  determined  by  photo- 
graph comparisons. 

The  problem  connected  with  photograph  viewing  for  "n"  determinations 
is  that  effect  of  meandering  cannot  generally  be  determined.     If  one  is 
using  the  channel  meander  length  on  slope  determinations  then  the  "n" 
value  found  is  probably  correct.    However,  if  the  flood  plain  length  is 
being  used  (which  is  the  case  for  overbank  flood  flows  )  "n"  must  be  ad- 
justed for  the  degree  of  meandering. 

For  the  in-channel  flow  condition,  "n"  was  determined  as  discussed 
above.    The  meander  channel  slope  was  used  for  this  condition.  However, 
for  the  overbank  flow  condition,  "n"  on  the  flood  plain  was  determined  by 
previously  explained  methods  and  was  higher  than  the  channel  "n"  used  for 
the  in-channel  condition.    The  slope  for  the  overbank  flow  is,  of  course, 
a  steeper  one.     "n"  was  taken  to  be  O.OO5  higher  for  the  main  channel 
portion  of  the  cross-section  than  for  the  overbank  flow  case.    This  is 


dons  to  account  for  tha  fact  that  much  oi  the  fiov/  in  the  area  of  the 
channel  Is  actually  fol,*loi;iDg  the  channalo    Ihic  .is  core  convenient  than 
assuming  that  ths  channel  portion  of  the  flov?  is  acting  on  a  jnilder  bIqj.q^ 
S5i£?i!.'--'^-ij^  Often  the  eroGS-soction  of  a  natural  stz^eara  is  so 

irregular  that  the  I-^nniDg  equation  cannot  be  conveniently  applied  to  the 
crosD-sect ion  as  a  whole.    For  ths  nster  surface- elevations  ^;hers  tb3  flcif 
is  in  transition  bet\';'ean  chann;^!  and  flood  plain  i'loj j  neither  the  wetted. 
I>erirn,eter  P  .nor  the  hydraulic  radius  E  scemratsly  reprccGnts  the  act^ial 
condition^    An  alternate  vjay  to  express  this  idea  io  to  say  that  the  ex- 
ponent on  E  clOuiR  not  ecnal  0^66"]  in  the  transition  zone  t,    This  is  eS" 
pecially  true  if  ths  rnin  chanr^^l  is  deep  and  tne  flood  plain  is  quite 
flat.    To  alleviate  this  problem,  the  cross ■■- s e ct ioiial  area  can  be  brohen 
up  into  several  portionn  roaking  the  portion  divisions  at  the  sharp  breaks 
in  the  section, 

CoiHputer  So.lution  of  M3.nnin5;'5  EquatJ.ons    A  conipntGr  progz-'am  ^jas 
written  by  ths  author  in  the  Fortran  II  language  for  the  EjM  3.620  com- 
puter (Appendix  I))  to  solve  Banning ' s  Equation  for  Q  at  a  given  croas- 
section  for  all  depths  (by  OiVi -foot  incrensnts)  betneen  specified  limits. 

Where  the  cross -sect ion  is  irregular,  ths  prograra  first  computes 
area  and  ^^'ottod  pariciater  for  each  portion  of  the  croGS-section«  The 
stage-discharge  information  is  then  pi'inted  out  for  that  portion*  The 
same  calculations  are  rnade  and  printed  out  for  all  the  portions  of  a 
particular  cross-sect ioue     Finally^  the  stage-discharge  values  for  all 
portio'is  are  sunsaed  up  and  printed  out  as  tlie  stage-discharge  relationship 


valid  for  the  ^^liole  cross-sactione     This  suisnary  informtion  is  then 
used  ar?  input,  to  the  PF?o 

ge stills  th@  Reliability  of  i^lBEHioIJ!  Sguaticns    The  reliability  of 
Manning's  equation  and  of  the  coniputer  program  was  tf-sted  by  coapai.'in^ 
the  st a/i;G -cl iscbarge  rating  curvos  'i^ith  thoys  \.hich  have  be^ia  develcpea 
by  the  U3G3  for  the  East  Forh  (Figuj-e  10)  of  Duck  Creek  and  for  Duck 
Creek  (Figure  11 ) ,    Tae  rating  curvG  for  the  East  Fork  was  dove loped  for 
the  saine  location  (previously  described)  used  in  this  study ^    The  rating 
curve  for  Duck  Creek  as  developed  by  'fiie  USGS  is  based  on  thair  staff 
gage  located  at  the  bridge  uhere  back  Creak  flous  under  the  county  road 
between  sections  11  and  Ik,  Tj-TH^  F^ioE.    The  l.rontana  State  University 
stage  recorder  was  originally  placed  at  this  sanie  Iccation  but  was  ivoved 
in  June ,  19 6 '4  to  a  point  330  feat  north  (down.strcam)  of  the  county  road 
centerlineo    Tlie  old  rating  curve  does  not  now  apply  to  the  new  location 
of  the  I'Bll  .stage  recori?er  and  since  this  paper  iucludes  a  study  of  the 
runoff  event  of  June  .I6,  I965  a  new  rating  curve  had  to  be  developed « 

Flow  has  been  meacured  on  three  cocas ioDS  by  the  USGS  since  the  re- 
corder on  "Duck  Creek  was  moved,    Tv^o  of  these  measurements  were  E:ade  in 
the  spring  when  the  flow  was  over  an  iced  channels    Therefore ^  an  acciu'ate 
rating  curve  cannot  be  drawn  based  on  actual  laeasureiirants »     It  was  de- 
terunnad  that  the  F^anning  equation  gave  a  surficiont3.y  accurate  rating 
curve «.    This  was  done  by  comparing  the  actual  and  calculated  rating  curves 
at  the  recoruer  on  the  East  Forkc    At  this  location_,  the  rating  curve 
given  by  Laaning'c  equation  fell  near3y  on  top  of  the  best  visually  fitted 


-J  1  ,  I  I  I  !  I 

LTv  j^-  OO  CVJ  HO 

Outside  Gage  Reading  (feet) 
Figure  10;    East  Fork  of  Duck  Creek  Rating  Curve 
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Outside  Gage  Reading  (feet) 
Figure  11:    Duck  Creek  Rating  Curve 


line  trirough  pointc  Lisaeured  by  the  USGS*    Since  this  is  true^  it.  may 

he  assUiiied  thau  tho  rating  curve  for  Luck  Creek  caieulated  by  the  I-feDning 

rulationyhip  should  "be  adequatcs  for  pi-rpcses  of  bhis  stu^jo 

As  already  stated^,  t"J0  different  slopes  have  bson  utilizod  for 
developing  rating  curve?  using  wannliig'G  equation,  depericling  on  vhcthi^r 
the  flov/  is  mainly  in-cbannal  or  on  tha  flood  plaiUo     In  the  forcisr  caso, 
the  length  of  the  n-ain  channel  center  line  or  .the  raander  lergth  was  used 
to  determine  the  bed  slope »     In  the  latter  casc^  houever,  the  length  of 
the  flood  plain  ^;ao  used.     This  length  is  of  course  shorter  than  the  meand 
length  giving  the  flaodplain  a  greater  slope  than  the  re::in  channels  In 
the  case  of  Rich  Creeh,  the  average  ratio  of  iiieander  to  flcodplain  length 
is  about  lc89c 

Curve  Numbers  and  Antecedent  |;;Od£ture_  Condition 

As  sho^m  in  Chapter  III^  the  soil-cover  complex  number  (Crl)  plays  a 
vital  role  in  the  SCS  method  of  detersd-ning  runoff  e    A  eoir^posite  CN  muGt 
be  calculated  for  a  ^"atershed  if  more  than  one  soil  type  or  cover  con- 
dition exist;; e    The  curve  numbers  uced  by  the  SCS  for  various  combinations 
of  soil  and  cover  conditions  are  given  in  Appendix  C« 

A  conpcoite  Qj-I  was  calculated  for  each  of  the  2h  subnatersheds  on  the 
Duck  Creek  watershed.    The  method  used  to  co-npute  this  co^nposite  is  shov^n 
beloi/  and  is  sumna-rized  in  Table  3=     First  it  Kas  determined  \?hat  the  CH 
v?ould  be  if  the  v/hole  subv;atershed  consisted  of  one  soil  type^    Tliis  nas 
done  on  a  weighted  basis  according  to  the  percentages  of  acres  falling  in 
the  various  cover  classes*    Tliis  was  done  for  each  soil  type  (typ.ss  B  and 


(1) 

Soil 

'  (2) 

Cover 
Condition 

(3) 

Curva  Number 
frc;m  AppeDclix  C 

(•'■0 

^  of 
total  area 

(5) 

Column 
(3)  X  (k) 

(6) 

%  of  area 
of  Bull  typQ 

(5)  X  (6) 

B 

Strip 
fanniDg 
Ran^'e 

80ei5-'^'              21, h  2:(o2 
CH  if  all  typs  B  soil  71.i:- 

29.3 

r-0.9  : 

 „j 

C 

Strip 
far-rij.i-g 
Eaiige 

67-                   21.k  18.6 
79                 78.6  62.0 

70/r 

57^0 

J 

CoriipOGite  CN  (roun<5.cd  to  nearest  yholo  nuru-osr )  » 

1 

78  ! 

! 

i 

Table  3  J    Calculation  of  Conposite  Curve  Ilircibsr 
on  a  Sam-ola  Subvjatersbsd 


•^n  average  of  CN's  for  fallov?  and  STiiall  grain  (straicbt  re;  in  good 
hydrologic  condition)  covur  condition  classes « 


C  vere  the  only  typas  occurring  on  Duck  Creak  )• 

The  v?eighted  CS's  oboal}i:?d  for  each  soil  type  ind5.¥idiial.ly  are 
then  again  ueightsd  v/lth  regard  to  the  percent-ages  of  area  i?ithi)i  each 
Biibwaterohed  that  fall  into  the  various  coil  groups »    ThiD  rQf:thod  of 
figuring  a  coi-posite  CN  is  based  on  the  assUiSption  that  the  various  soil 
type  porcentages  used  occur  in  those  cain.3  percentages  on  each  of  various 
cover  condition  groups  used-    If  this  ia  not  true ,  L-ign  if  leant  errors  in 
Oil  calculation  could  occur »     It  is  felt  that  such  an  error  uould  not  be 
irfjportant  in  the  cas«  of  Duclc  Creek  since  in  noarly  all  Eub\/aterGhodSj, 
one  f-oil  type  or  cover  condition  v?as  prsdorainant «    Wnen  this  is  not  trv.-3 , 
it  is  nscecsary  to  calculate  the  v?aighted  CN  for  each  subarea  of  each 
Eubivatsrshed ,  i:hereupon  the  .breakdown  of  r>oiJ.  type 3  for  each  su.barea  vould 
have  to  be  dcterLained,    This  detail  waB  not  considered  necessary  on  Duck 
Creek* 

The  determination  of  CN  from  Appendix  C,  ^vhich  is  based,  on  an  average 
antecedent  u.oisture  condition  {Al.:c)_5        input  to  the  PFP^    Also,  the 
actual  AI'^C  is  deternined,  as  described  in  Chapter  111^  and  read  inbo  the 
PFP«    Tno  adju^jtmenb  of  CN  is  internal  to  the  r??»    iOr  the  stom  of  June 
16,  3.965  on  Duck  Creek,  the  AI^C  vas  determined  to  be  dry  ('Cype  I), 
Tirie  of  Conccrrfcrati^n  • 

Tae  time  of  concentration  for  an^'  given  cheDnel  reach  is  taken  to 
be  the  tirr.e  of  travel  for  full  bank  flow  conditions.  The  velocity  for 
full  bank  flow  is  doteru?ined  by  Manning's  equation  at  both  ends  of  the 
channel  reach.    Travel  time  is  then  given  by  the  follouiug  relxitionship : 
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o  r 

L  is  the  reach  leDgtli  in  feet,  V-j  ard  Vp  arc  the  velocities  at  the  in- 
flow and  outflow  sections  of  the  channel  reach  and  t  is  the  travel  time 
through  the  reach c     In  the  case  uhors  the  reach  in  quest  ion  11;  taetv^ssn 
a  etreaiii  cross  section  and  the  boundaiy  of  the  watershed,  velocity  is 
taken  to  be  zero  at  ths  boundary.    These  times  are  then  added  down  tbroi'3h 
the  watershed  to  find  the  vzatershed  lag  tine  (tiine  between  centers  of  mass 
of  the  rainstorm  and  the  outflov?  hydrograph )  for  the  waters  he  do    When  this 
calculation  does  not  match  the  acbual  lag  ti^ui;  as  observed  frcni  outflow 
hydrographs,  each  reach  travel  time  must  be  adjusted  by  proportional 
multiplication  so  that  the  actual  lag  tiine  is  obtained.    The  tiria  of  con- 
centration for  each  individual  sub'?atershed  is  then  taken  to  be  a  portion 
of  the  stream  full-bank  flow  time  plus  overland  and  tributary  channel 
flow  times The  conibination  of  laain  stream  and  side  channel  travel  tiries 
that  is  longest,  is  taken  to  be  the  time  of  concentration  for  the  sub- 
watershed* 

The  tributary  channel  travel  times  were  based  on  the  velocity  of 
two  or  three  feet  per  second  at  the  mouth  and  zero  at  the  head^  These 
two  velocities  and  the  length  of  the  tributary  channel  are  then  sufficient 
to  give  the  tine  of  travel  as  given  by  equation  (35 )«    OTlrls  is  an  approxi- 
mate method  but  sin^T iciently  accux^ate  on  the  Duck  Creek  subwatercheds  be- 
cause they  are  characteristically  long  and  narrow.    The  tributary  channel 
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timec  of  travel  uore  short  compared  to  the  mia  channel  travel  times ^ 
For  shorter,  \-nd.er  sub'viatersheds  the  determination  of  time  of  t  rave  J. 
and  time  of  concentration  would  noed  more  scrutiny,  for  the  side  channel 
tlmoE  would  const itutfe  larger  percentages  of  the  total  timaSc . 
Storm  Rj^ta 

Both  actual  and  design  storr-is  can  be  handled  by  the  Pfp,    The  rain- 
fall amcunt  is  always  read  in  as  either  the  actual  aniount  froiii  scire  storni 
or  the  design  rainfall  s-inount  as  dictated  by  scire  design  criteria^  The 
pattern  of  an  actual  stoma  can  be  taken  into  account o    For  instance,  whan 
the  watershed  is  broken  down  into  several  subwatersheds  a  different  rain- 
fall amount  can  be  designated  for  each  watershed „    Temporal  variations  can 
also  be  accounbed  for  in  actual  Etorj;:S,     The  starting  timi;  of  the  rain- 
fall can  be  different  for  each  subwatershedc 

Rainfall  patterns  can  be  read  into  the  program  as  cumulative  dlmen- 
sionless  graphs  or  gx'aphs  with  the  actual  time  and  rainfall  amounts  as 
abscissa  and  ordinate,  respect Ivclye     As  many  as  nine  different  patterns 
can  be  used  in  any  one  run  of  the  computer  program, 

Tlie  amount  of  rainL'all  falling  on  each  of  the  2k  swhwatersheds  on 
Duck  Creek  was  determined  by  the  Thiesssn  ^lethod.    The  Thisssen  polygon 
was  drawn  for  the  seven  rain  gages  on  or  near  the  watersiied.    VJith  this 
method  the  rainfall  over  the  entire  Thlessen  subarca  is  considered  to  be 
uniform  and  equal  in  value  to  the  amount  recorded  at  the  gage  which  gen- 
erally falls  in  the  center  of  the  polygon  if  the  gages  are  equally  spaced c 
As  would  be  expected,  some  of  the  subwatersheds  fall  in  two  or  thi'ee  of 
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the  ThieGsen  'golyiionn ^  in  whi.ch  case  the  raiDiall  ariount  is  a  v/eightcd 
amount  based  on  the  pei-ce'ntages  of .  area  falling  in  each  of  the  polygons » 

Fortunately,  ilie  edges  of  tho  polygons  fall  close  to  the  ridge  bo- 
tusen  the  East  and  West  Forks  of  Duck  Croeh;  this  eliminates  the  largest 
criticism  of  the  Thiessen  method  vhica  says  that  the  method  does  not  take 
into  account  geographic  features ^ 

Storm  moveiiient  v/as  taken  into  account}  by  varying  the  starting  time 
for  each  of  the  snbv?atersheds ,    !I;he  speed  of  the  etorn  was  determined  by 
examination  of  the  records  of  the  three  continuous  recording  rain  gages . 
The  time  for  inOYeVLsnt  of  the  storm  across  the  -Yatershed  uas  found  to  be 
about  3-13  hours  for  tha  sstorm  of  June  l6,  j.965»    Tne  starting  tiineo  for 
the  various  eub^atersheds  I'ere  then  determined  to  the  neare£;t  0^1  hour  by 
visually  determining  the  relative  poco.tiona  of  their  ce/jtroidSe 
Sumaiary 

The  amount  of  and  kinds  of  data  supplied  to  the  PFP  vary  from  v;ater- 
shed  to  v?atershedc    Not  all  types  of  data  that  could  be  processed  by  the 
program  need  be  used  in  every  uatershed  model  developsd  for  the  coniputer 
solution.     For  example,  although  there  are  a  nui'iber  of  small  stock  pond 
reservoirs  on  Luck  Creek  uatershed,  their  effect  is  felt  to  be  minor,  and 
therefore  the  KOdel  does  not  include  any  reservoir  inforinaticnc     The  data 
supplied  to  PFP  has  been  discussed  in  this  chapter .    With  this  infornation, 
the  PFP  is  capable  of  calculating  the  runoff  hydrograph  frorn  Duck  Creek 
based  on  the  usual  SCS  hydrologic  formulas  and  procedures^ 


Chapter  yi 
ACTUAL  AlID  IIYPOTHSTICAL  STOIilB 
The  actual  otorm  of  June         I965  on  Duck  Croek  v/as  used  to 
ad^iiot  a  rainfali-runcjff  ticdel  of  the  vater^liiid,    G?£:en  hypothetical 
Btoiu-a  based  on  U„        V'eathor  Bureau  frequency  data  were  routed  through 
the  ricclal  to  predict  poak  discharges  <> 
Actuaj.  Stoi-iTi  of  Jurie  2.6 ^  ^Sp'^- 

:  ^  /   The  actual  storn*  used  for  model  adjustiisnnt  on  Duck  Creek  hat- 
he rat of or Q  been  called  the  stor^a  of  duns  lb,  1955,    Tbs  resulting  psak 
discharge  actually  occurred  on  this  date,  although  the  sr'corw.  ia  question 
start^^d  at  about  6  a<,n,  on  Jur-c  13  and  lasted  until  about  k  a^iiu^  duija  I6 
for  a  total  duration  of  JO  hours «    These  ti~nc?s  are  average  times  as  the 
stora  had  G-lightly  varying  durations  as  it  moved  across  ths  uatershede 

The  iiajority  of  rainfall  occurred  during  three  rain  burLits  approxi- 
matcuy  a  day  apart; «  The  duration  of  the  bursts  uere  approximat e ly  3;  ^-^ 
and  6  hoiu'a » 

Figure  22  shous  tlie  areal  dictribirtion  of  the  total  ariourt  of  rain-= 
fall  in  the  foi'm  of  an  isohyetal  mpo  This  is  based  on  the  records  frora 
a  total  of  seven  rain  gages  on  or  near  the  ijatershede 

The  stora  liioved  generally  in  a  northrestsrly  direction  dcun  the 
watershed's  main  channal.    An  indication  of  the  stora  is^ovc-aanb  vas  ob~ 
tain^'d  by  studying  records  of  time  versus  cimulativ-a  rainfall  amount  by 
percentage  for  tuo  of  the  three  contiDuous  reeordirg  rain  gagr^s^  These 
indicated  that  the  difference  in  tiine  at  the  tuo  recorder  points  at  yhich 
the  soiiie  perce^rbages  of  total  storn  amount  had  occurred  v/as  about  tuo 
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a  tfep  ox  Juri^i  16,  1965  St03:iii  on  Duck  Creak  V/a';eriJhc;d 


hours  6    This  was  approziuiately  trus  for  all  three  bursts  so  that  on^ 
storm  pattern  was  assumed  to  bo  adequate  for  all  of  the  si>>)':Jatar:^liedG  e 
Assumiiig  that  the  storai  traveled  at  a  constaril}  rata  of  spsed.^  tha  2s0 
Yiom:  difference  in  tiina  of  storm  beginning  "betueen  the  contizi-uiou's  re- 
cordiris  points  becomes  3^5  iiours  from  the  areal  cerrbrold.s  of  the  tv?G 
K03t  yxd2ly  separated  siibuatarshads  o 
JilX)t  haMca  1  Stor^ 

Storm  durations  of  12  an<3.  2h  Lours  viere  cho:'en  arbitrarily  for 
the  purpose  of  using  the  5C3  i:ii:thol  to  predict  poak  di-chargeSo  BaiD- 
fall  amo rants  routed  for  both  durations  were  in  0,5  inch  increriie-nts  Pnd 
over  a  range  si:ch  that  5,  10^  25/  50;  ai'.d  100-year  frequency  storms  would 
ba  included 0    These  amounts  are  thore  predicted  by  ths  VJeather  Bureau  in 
their  Technical  Paper  Mo*  ^-i-O  aa  prepared  by  favid  Ilersbfield  (1951). 
Therefore^  for  12-hour  ctoraS;  pracipitat ion  P  of  2e0^  2^5^  3.0^  3c 5^ 
and  hcO  inchsGi  and  2^i-"hour  storrjs  v?ith  a  P  of  2,0^  2.5^  3<.0;  3o5^  kcO^ 
and  k^y  inchas  v'ore  rout'^dc, 

The  time  of  storm  beginning  was  taken  to  be  constant  on  all  cub- 
v?atarKh3ds  sinc2  there  is  little  iiif orj^iation  on  the  direction  of  pre- 
doiBinant  stoTm  riOvremant  on  Dack  Creek c     If  a  coniiistont  pattern  of  storn 
movement  \?as  kno'^^n  in  the  case  of  a  particular  watershed,  then  some  varied 
time  of  storra  beginning  could  be  aBsuraed;  even  for  h;^^othetical  BtorriSs 

Also_j  the  total  amount  of  rainfall  imst  be  assrred  to  ba  constairfc 
on  all  sui^watershsds  unless  there  is  inf orEiation  irdxcating  that  the 
Btorm  center  for  most  knox-m  actual  storms  Is  consistently  located  in  one 
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portion  of  the  watsrshcdo    Suca  Information  is  not  knoun  for  Duck  Croek« 
Therefore?;,  rainfall  aY^iomits  uero  conGldGred  to  be  spa^tlallj  InvarirMjl  „ 
The  effect  of  tiiL  location  cf  storm  centers  could  be  ctudied^  hov?Qver=, 
hy  cantering  storms  in  varS.cii^  parts  of  the  ijatcrohed  and  coDipariDg  the 
calcu3.at$;l  o^itflov^  hyclrograpb.  results »  ■  . 

The  actual  runoff  h^'^drosi'aph  for  the  storra  of  June  l6.  1965  is 
shov?n  In  Fxgura  13^    TIiIb  hydrograph  had  baen  devslop^;!  by  previous  in- 
vestigatoris  bascd  on  the  best  dateraiiiation  of  the  rating  curve  iu  effect 
at  the  nev;  ctroam  otage  rscoruer  locaticno    The  voliiss  of  runoff  associ- 
ated uith  this  cur\^c  is  0.^.3  inches  for  the        square  milo  Duck  Creek 
vjatershedo 

After  thic  study  \;&3  i/ell  under  '-Jay  it  v?as  discovered  that  tuo 
errors  had  been  cojiomitted  in  the  reduction  of  the  continuous  sta^e  ra« 
c order  data.    Tha  first,  made  by  studvnts  v? or king  on  the  Drainage  Corre- 
lation Research  Project,  was  that  the  vol'.v-2?.  previously  calculated  vac- 
in  error e    Tha  resulting  yield  should  have  been  0.23  instead  of  the  0.^3 
inches c    The  second  error  was  r,:ade  by  the  U.  So  Geological  Survey  in 
their  publications ,    The  peak  stage  printed  for  the  Jure  l6^  19^5  J'^torm 
van  k^'jh  feet  v?bich  gives  a  peak  of  87 0  cfs  based  on  their  rating  curve 
for  Duck  Creek o     The  stage  value  should  have  been  ^*.68  feat  uhich  gives  a 
peak  of  jhO  cfs  instead.     Since  the  original  rating  curve  for  the  njn? 
locatioij  v?as  based  primarily  on  this  one  storm,  it  is  nov?  invalids 
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The  author  lias  clcvelopsd-*^  another  ratiD^  ci-rve  based  on  the  tvSMining 
Equation^  this  bsi^g  justified  for  reasom  described  in  Chapter  Yo  Onlj 
two  points  are  kno^;n  on  the  rating  curve ^     Oae  Ib  the  stcria  being  tested 
and  the  other       the  elevation  or  stage  at  i^bich  the  floi?  is  i::ero»  A 
rating  curve  uas  dra\-ni  throu{j5b  these  ti?o  point:;  based  on  tl:s  J/aiining 
EquatiGiio    The  runoff  bycirograph  ^vas  again  draun  based  on  the  iisvi'  ratlDg 
curve ;  anQ  is  a3.G0  slic^vvn  in  Fi^jUre  I3*    The  ordinates  of  this  plot  vers 
read  into  a  computer  program  \-^ritten  by  Gary  I3^7i3  to  deteririina  the  yields 
Tlie  yield  nor?  believed  correct  for  the  June  I6,  I965  storm  on  Duck  Creek 
is  Oc3l^  incbeCs 
gopip liter  I3vm£  Com]?leted 

A  nuraljer  of  runs  or  storra  routings  v?ars  made  on  t^o  different  ccca- 
B ions .    These  runs  are  sumirarised  in  Table  4. 

Runs  1-A  through  1-M  and  2~A  through  2-M  are  identical  except  for 
th3  curve  nuuibera  (CH)  uocdc    The  number  1  runs  make  use  of  a  set  of 
en's  originally  thought  to  bo  those  in  effect  on  Duck  Creek o    For  the 
number  2  runs,  each  CI'J  of  the  set  vac  reduced  by  the  value  of  three » 

The  A  ana  B  runs  (both  nuabers  1  and  2)  caused  the  June  I6,  1965 
storm  to  be  routed  through  the  Dack  Creek  modele    The  only  difference 
bctvje'sn  A  and  B  is  the  assumption  icade  about  the  time  of  the  storm's 
begiDuii^go      ■  '   '  '  .       .  ■ 

Runs  C  through  lA  are  for  hypothetical  storir.Sc    Kuns  C  to  H  are 

■^With  assistance  from  Gary  lev  is ^  graduate  student  at  Montana  State 
Unive-:-\sity  uorking  t0;-/ards  an  M,  S«  Degree  in  Civil  Engineering. 
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fabln  k :    Suri~ary  of  Coriputer  Runs  Completed « 


2'i-b.our  duration  ptcnns  i/ith  variable  rainfalJ.  anovmis ^  tjhile  I  to  H 
are  12-haur  dui'ation  stor^is  yitli  variable  rainfall  EiuountSo 


■    ■  Chapter  VII 

RESUKTS 

Tiie  data  available  on  'the  stora  of  June  'X6^  l?^^?  on  the  Duck 
Creek  watershed  was  used  to  set  up  the  computer  model  for  analysis  by 
tba  Soil  Conservation  Service      Project  Formulation  Program  (SC3-PFP}« 
Tv?o  sets  of  coiiputer  runs  Mere  made  by  the  SC3  in  liyattiiVilla I^^iryland 
on  an  IBM  ^oO/hO  ole^ctronic  cOiitjuter. 
Actial  Storm  SvajAjayon 

Tiie  initial  rim  gave  a  total  watershed  yield  of  0^57  incheSo  Ojhis 
v.* as  higher  than  the  actual  yield  measured  at  tho  sx-reaTn  gage  which  was 
0,35  inches.     However^  at  the  t ime  the  eonputer  rur-s  v/ere  Eade  the  actu 
yield  vas  thou{rht  to  ot  inches'  (due  to  previously  explairad  error) 

The  model  yas  initially  set  up  by  evaluating  the  soil-cover  com- 
plex nuaibsrs  {CN)  thou^cht  to  be  thooe  in  effect  on  Puclc  Creek «  These 
are  based  on  the  chart  used  by  the  SCS  (see  Appendix  C  of  this  paper). 
The  procedure  used  by  the  SCS  is  to  check  the  yield  for  the  originally 
determined  set  of  CN's,    The  Cil's  are  than  adjusted  until  t]ie  calculate 
yield  equals  the  actual  yield.    This  is  considered  to  be  the  calculated 
runoff  hydrograph  which  can  ncj  be  compared  to  an  actual  runoff  hydro- 
graph  if  it  is  available e  ■ 

In  an  effort  to  correct  the  yield  to  Oji3  inches  (thought  at  the 
time  to  have  been  the  actual  yield)  the  author  reduced  the  CK's  for 
all  sul;jnater:;heds  by  a  value  of  three  and  the  second  set  of  couipncer 
runs  was  r/.ade .    The  resulting  output  uas  hydrographs  having  yields  ex- 
actly equal  to  OA-'S  inches « 


CompirfcGr  Hesults  i    Figure  1^^  shoi/s  the  ac tua],  runoff  hydro^;};aph 
for  the  June  l6,  lyS^  storm ,  points  indicating  t'pe  computed  hydrograpb. 
peo.Us^  and  t\JO  liydrojraphs  vjhxch  have  baen  postulated  on  tha  basic  of 
the  best  evidence  available c    The  actual  hydrograph  as  measured  has 
a  lov7er_^  J^onger  profile  than  the  calculated  hydrographc;  ^    The  cal- 
culated hydrograplia  Lave  higli  peakG  ard  stcepar  recassion  limbs o 

The  f oin;'  inolatod  points  shoyn  on  Figure  1^  are  the  highest  peak 
discharges  calculated  by  tha'  computer  runs  1~A^  l-'B^  and  The 

first  runS;,  each  with  a  yield  of  0«57  ixichos,  are  indicated  by  the  tuo 
highest  points  (i-A  and  l~Ij)o    The  hi^pjer  of  these  (l-A)  io  for  the  case 
vhere  the  ticie  of  storm  beginning  va.s  assuKsd  to  vary  from  one  suo- 
watershed  to  the  next «    Tair-i  hydrograph  has  a  peak  discharge  of  I'l-lo 
cfs»    The  lover  of  the  two  (l-B)  with  a  discharge  of  22.h^  cfs  is  based 
on  the  assufnption  that  the  starting  tii^a  \ms  the  saraa  for  all  siib;jatcr- 
sheds c    Tile  later  comp^rter  runs  x^ith  reduced  curve  nuiri^^fcr  values  re- 
duced tha  yield  to  OJ13  inches ;  and  the  peaks  (2"A  and  r!-i3)  to  llhh  and 
1006  cfs  resp2ct iva3,ye    Agaixi  the  higher  peak  vas  for  the  case  of  the 
variable  storra  beginnings 

The  tifo  hydrographs  shoun  in  th^^ir  coiaplets  form  on  Figure  1^  have 
been  postulated  by  the  author,  and  are  based  on  the  results  of  the  last 
coinputer  rurs^     Each  of  these  hydrographs  has  a  yield  of  Oc35  inchss, 
matching  the  yield  of  the  actual  hydrography    Their  shapes  ore  identical 
to  those  of  the  later  ccii^putcr  rune*     If  further  computer  runs  should  be 
made  with  loivar  Cli  values,  it  is  expected  that  the  shapes  of  the  hydro- 
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graphs  DiigiVij  chr'Dge' s light ly»    Coiriparlson  of  th^^  coiiiplete  "hydrograx^hs 
of .  the  fom'  computer  runs  tvhlch  vjere  rrsde  indicates  that  the  change  In 
shape  t-?ould  not  b-s  great  _j  hoy  ever* 

As  can  dq  seen,  there  is  a  large  discrepancy  bet-jeeii  the  calculated 
and  act\ial  hydrographs  with  re^-ard  to  peak.     Ths  SC3  states  that  tha 
method  was  Dot  designed  to  ma-Gch  the  Iiydrcgraph  in  shape  but  only  in 
paak  value «     Hovever,  in  this  casa,  the  calculated  pual:  is  either  2^kl 
or  2cT4  times  as  great  as  tha  actual  peak  depending  on  vhich  assumption 
is  used  for  the  tins  of  siors  beginnirig  on  each  sub-ivatei'shed^    The  rami- 
fications of  these  results  will  be  discussed  in  Chapter  VIXI« 
Hypolnietical  Stcrra  Evalp.atlon 

Ona  purpose  of  this  study  uas  to  precUct  flood  peaks  for  various 
frequencies o     Ihis  \?as  accomplished  by  routing  a  number  of  hypobhstical 
storijis  through  the  v/atershsd  model  to  determina  the  paak  discharge  pre- 
dicted by  the  SOS  risthod.    Kie  final  results  of  this  portion  of  the  work 
Is  suriiinarized  In  Figure  I5,    Of  course  no  actual  stornis  have  occurred  by 
which  these  discharges  can  be  compared. 

Figure  I5  is  tertative  and  should  not  be  used  for  design  purposes e 
It  is  shov?n  here  only  as  a  deraonstration  of  how  the  SC3-PFP  results  of 
routing  of  hypothetical  storiiis  can  be  used  in  conjunction  with  the  S» 
Weather  Buresu  frequsncy  data  to  predict  peak  discharges  from  S]rBll  v'ater 
sheds.    The  two  curves  narked  \JB  (for  Weather  Bureau)  are  plots  of  ex- 
pected precipitation  P  versus  the  return  period  in  years  (abscissa  scale 
at  the  bottom  of  the  graph)  for  the  storn  of  magnitude  P  for  both  2h-- 
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Figure  15:     Example  of  Peak  .Discharge  Design  Chart 


and  X2-hour  storviis  as  indicated.    The  other  tuo  curves  r:arV;ecl  DC  (for 
Duck  Cre3k.)  sho:;  th.3  p-Xcts^  of  F  versus  paal:  dxfichnrga  q-Q  (abscissa 
scale  at  the  top  of  th^^  graph)  for  the  Di;.ek  Creeh  uatcr&hsd  as  deter- 
mined by  the  computer  nod  el  us  1113  SCB  h5'"d:':'ol03lc  dasign  critGria  for 
the  Game  t^'JO  stem  durationSc     These  ttvo  curves  represent  the  Liodel  as 
set  up  for  tha  later  coaputcr  runs  (not  as  yat  adjusted  as  earlier  e3^■■ 
plained  and  thuc  not  suitable  for  design  use ) «    If  add it i oral  adjustment 
of  the  tiodel  could  "be  imde^  the  same  procadtire  as  descrihed  here  would 
he  applicable  in  setting  up  a  siailar  set  of  curves  for  use  in  dssigUc 

If  Figure  I5  Mere  to  be  used  for  d£sign_j  the  follo^ving  i^roQii&DXQ 
should  he  follo!f2d  to  properly  predict  the  p^^ah  dischar{^e  for  a  ftorvAc 
Generally,  structures  on  a  uatersaQd  are  hydraulically  dasigned  to  func- 
tion  properly  for  some  peak  discharge  cau3Gd  by  a  storr?!  of  a  particular 
frequency  or  retiurn  pariode     Suppose  for  an  exarple;^  as  f'houn  on  Figure 
15j        plan  to  deaign  a  culvert  that  will  handle  a  ^O-ycar  frequency 
storni  of  a  duration  of  2h  hours.    The  procedm^e  is  to  start  at  the  bottoci 
of  Figure  I5  v;ith  the  ^0--yeP.T  storm  and  proceed  verfcically  along  the  lira 
sho'/n  until  the  hour  curve  is  inato    Then  if  the  ordinate  P  is  read, 

ue  find  that  the  Weather  Bureau  predicts  that  the  precipitation  is  3,75 
inches  for  a  50»;year  -  2^^-hour  dui'ation  etoriiic     If  honever,  we  follow  a 
horizontal  line  (to  the  right  in  this  example  )  liutil  we  meet  the  dis- 
charge for  a  2^:™hour  storm  for  JJ.i.oh  Creek;  and  then  proceed  vertically^ 
we  find  that  qp  ^  320O  cffa. 

Ga  the  other  hand,  if  the  peak  discharge  for  a  runoff  event  had  been 


„  77  - 

measured  on  Duck  Creek;  the  exact  opposite  procedure  to  that  dos- 
cribed  in  the  previous  paragrapbj  could  be  folloued  to  determine  the 
frequancy  or  reti^rn  psriod  of  that  stormc 


.  Chapter  VIII 
Discussion 

The  investigation  described  in  the  preceding  chapters  has  been 
of  two  types.     First a  model  of  Duck  Creek  watershed  vjas  developed 
and  then  tested  by  using,  as  input  data  to  the  computer  program,  in- 
formation describing  an  actual  storm  which  occurred  on  the  watershed 
on  June  l6,  I965.    The  resulting  synthetic  hydrograph  generated  by  the 
program  was  compared  with  the  actual  hydrograph  recorded  at  the  mouth  of 
the  watershed.     The  model  was  adjusted,  after  initial  comijarison,  and  the 
program  re-run.    Second,  the  Dack  Creek  model  was  used  to  predict  the 
hydrograph  from  certain  design-type  storms  based  on  Weather  Bureau  storm 
frequency  data. 
The  Storm  of  June  I6,  I965 

The  results  from  the  computer  runs  for  the  June  I6,  I965  storm  do 
not  compare  very  closely  with  the  actual  hydrograph  recorded  for  the 
event.    The  discrepancy  is  thought  to  be  mainly  the  result  of  imper- 
fections in  the  model,  but  there  is,  however,  a  possibility  that  the 
actual  hydrograph  is  not  entirely  correct.    The  peak  discharge  that 
actually  occurred  is  felt  to  be  accurately  known,  because  it  was  de- 
termined from  the  USGS  rating  curve  for  their  crest-stage  gage.  This 
rating  curve  is  based  on  ten  current  meter  measurements  over  a  period 
of  at  least  ten  years.     Except  for  the  peak  discharge,  the  rest  of  the 
actual  hydrograph  was  determined  by  relating  depths  recorded  at  the 
water  level  recorder  to  discharges.    As  has  been  mentioned  earlier,  the 
water  level  recorder  is  located  several  hundred  feet  downstream  on  Duck 
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Creek  from  the  crsst -stage  gage,  and  the  stage -discharge  relationships 
for  the  two  stations  are  quite  different.    The  best  rating  curve  avail- 
able at  the  water  level  recorder  uses  ^fe.nning's  equation  and  is  drawn 
through  only  two  known  points. 

Possible  imperfections  in  the  model  seem  to  fall  in  one  of  the 
following  categories:    1)  preparation  of  the  input  data  and  2)  the  SCS 
method. 

Preparation  of  the  Input  Data ;    When  one  is  dealing  with  hy- 
drologic  problems  the  situation  never  occui's  that  there  is  an  over 
sufficiency  of  data  available.    There  is  always  a  limitation  on  the 
time -dependent  data  available,  and  an  economic  limitation  on  the  thor- 
oughness with  which  the  watershed  parameters  can  be  obtained.    A  rather 
significant  amount  of  watershed  data  were  obtained  for  the  study  re- 
ported herein,  including  the  results  of  a  medium  intensity  soil  survey, 
numerous  infiltration  tests,  a  complete  stadia-level  sujrvey  of  the 
principal  streams,  and  more  than  20  channel  cross -sect ion  surveys.  It 
would  seem  that  the  watershed  parameters  were  delineated,  therefore, 
more  completely  than  could  be  expected  in  the  typical  smaj.l  watershed 
design  investigation.    Results  of  the  study  show  that  even  this  pre- 
cision may  not  be  adequate  to  insure  great  accuracy  in  duplication  of 
an  actual  hydrograph. 

The  way  in  which  various  watershed  characteristics  were  utilized 
will  be  discussed  below. 
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The  curve  m^iibers  (CN)  vhich  are  v3gO.  in  the  SC3  program  to  charac- 
terise coil  uyp2  mud  land  use  treatment  \^er3  cieve loped  frora  large  qiian- 
tities  of  data  vhicb  were  obtained  from  all  parts  of  the  United  States, 
Sii?ce  they  represent  average  conditions,  and  ijsre  not  developed  £;peci- 
fically  for  I'ontana ,  it  coeiiiS  likely  that  u<50  of  the  CN  table  in  Appendix 
C  K;ay  be  responsible  for  some  of  the  discrepancy  in  the  reaultSo  In 
particular,  it  seems  that  the  adjn3tm2nt  nbich  is  made  for  dry  antecedent 
moisture  conditions  (type  l)  irsey  not  be  large  enough.    It  may  be  that 
for  Montana  watersheds,  the  adjustment  between  average  and  dry  conditions 
should  be  larger  thB;j  that  used  nat icnwide 

The  CN's  used  in  the  first  set  of  computer  runs  ^-u^re  adjusted  for 
the  second  set  of  runs  by  reducing  each  nuxiber  by  a  value  of  three .  Ihis 
reduction  theoretically  should  have  been  proportional  to  the  original  value 
of  each  CK  for  each  subvatershed.     It  is  uncertain  what  the  effect  of  re- 
ducing in  equal  instead  of  proportional  amounts  might  be. 

There  are  indications  that  the  tiraing  of  the  runoff  in  the  model 
vjas  in  error.    This  seens  to  be  especially  trii?  in  the  case  of  the  high- 
est peak  shc:-n  on  Figure  Ik,    As  can  be  seen  from  the  main  portion  of  the 
actra]  hydrograph,  lesser  peaks  occurred  both  before  and  after  the  rnaxi- 
laua  discharge „    The  calculated  hydrographs  seem  to  have  grouped  all  three 
peaks  into  one  and  corsequently  gave  too  high  an  estimate  for  the  mexi- 
mum  dischargee    "j^he  factor  that  controls  the  timing  on  the  model  is  the 
time  of  concentration  which  is  determined  for  each  sub-'atcrshed  =  As 
mentioned  earlier^  only  one  tiPie  of  concentration  csn  be  supplied  for 


-  81  ^. 

each  s ubwatersheo. ^  anfi.  tharafore  no  variatio;'.  can  be  aLiSU/iicd  in  this 
tiK3  for  different  discbar£^es ,    Tae  t  Ikg  of  concentration  can  be  sho'.;n 
to  vary  with  rate  of  clischarge  so  that  this  bocoHsas  a  limiting  factor 
in  the  construction  of  tho  nocielo     lha  times  of  conconbration  for  this 
study  v?ere  based  on  valocitics  as  calculated  by  the  ifenxiing  relationship 
for  full  bank  flovu    Tae  actual  event  iD  I96?;  prciucad  little  or  no 
flooding;  it  f:eeF.s  probable  therefore  thst  the  chan-aols  were  less  than 
bankfull  rnoKt  of  the  tine^    The  true  velocrc-ies  v/era  probably  las&  than 
'those  calculated,  anl  the  times  of  concentration  larger  than  thoso  Ccl- 
culatedo    This  eoens  lihaly  to  account  for  the  fact  that  the  calculated 
peaks  occur  at  tines  earlier  than  the  observed  peaks c 

The^  SCo  hsthod:    The  SCo  iiipthod  was  developed  before  eJ.ectroDic  corii- 
putei's  could  be  utilised  to  solve  such  problercs  and  therefore  had  to  be 
simply  constructed  so  that  it  could  be  solved  by  haii-d.    And  as  is  true  of 
a 3,1  methods  developed  for  predicting  peak  discharges,  the  SCS  oethod  is 
based  on  s^any  assir-ipt ions »    Assusiptions  are  rede  either  because  it  is 
thouglit  that  they  utll  not  greatly  affect  the  solution  or  because  they 
must  be  isade  to  irnke  the  solution  econoraically  possible  e     I^s^ny  times  in 
solving  hydrologic  problems ,  assumptions  niust  be  made  because  there  is  a 
scarcity  of  inf orr-ation  available,  even  though  it  is  recognized  that  these 
assumptions  my  significantly  affect  the  so3.ution« 

The  SC3  method  can  be  thougl^t  of  as  a  matheniat ica  1-graphic  model 
which  consists  of  a  number  of  formulae  and  graphs  o    A  .It  hough  some  flexi- 
bility has  been  built  inbo  the  ope?-^ation  of  the  PFP^  it  appears  that  the 


„  8?  - 

SCS  method  has  no'G  been  basically  altered  in  the  conversion  into  a  com- 
puter liiodel.  The  question  that  then  arises  is  whether  the  has io  •  nethoci 
itself  could  no\i  be  x£j^.de  nore  flexible  for  ccmpiri-er  so3..utlon3o 

Points  in  the  modal  that  cea:n  to  he  too  rigid  are  i    1)  initial 
abstraction,  2)  Icq  tir.if ^  3)  synthetic  unit  hydro^^raph  shape,  '■!■ )  storn 
pattern^,  &.nd  5)  '^he  reservoir  and  flood  routing  sathodSo 

The  initial  abstraction  L,  is  assumed  to  he  equal  to  Cn2  S  and  the 
lag  tiina  to  ba  equal  to  0.6  of  the  tims  of  concent  rat  ion«    The  coaffi- 
cients  of  Oc^^  and  0^6  are  certainly  for  the  avsragG  candid ionn  for  a 
mnnber  of  natersheds  that  vjcre  studied  by  the  SCS<,    The  rau  data  from 
which  these  coefficients  uere  determined  are  not  {^ivan  in  the  SOS's 
National  Sngineering  Handbook  (NEH),     Xf  these  data  had  been  supplied,  it 
might  be  possibj.e  to  choose  more  nearly  correct  factors  for  a  specific 
area* 

Also,  the  KSH  defines  only  one  synthetic  hydrograph  shape,  which  is 
intended  to  be  used  in  all  cases »    The  PFP,  ho'jever,  is  capable  of  absorb- 
ing any  unit  hydrograph  that  the  user  wishes  to  supply.    From  this,  it 
would  seem  that  the  SCS  method  as  outlined  in  the  lIEil  is  more  restrictive 
than  is  necessary  now  that  the  more  flexible  computer  programs  are  avail- 
able 6 

Only  two  possible  storm  patterns  are  given  by  the  SCS  for  use  in 
testing  hypothetical  storiiis^     Tij.e  pattern  most  connionly  used  (shonn  in 
Figure  $)  is  one  in  which  the  main  burst  of  the  storm  comes  at  about  the 
'midpoint  of  the  storm  in  tiiiis.    The  other  has  the  nain  storm  burst  come 


eailier  in  th":  storm,,    Again,  the  FBP_j  being  able  to  c?ccoiTiD.odate  aiiy 
etom  pattern  shs-pe ^  is  more  flexible  than  tba  Eiethcd  described  by  th^. 
SC3„ 

The  KSH  BU^(j;,eBts  sovoral  methods  of  fj.ood  routing  vjhicb  niay  ba  used^. 
Only  t'jo  of  these,  ho'jever,  have  be  on  incorporated  into  tha  computer  pro- 
gram^    The  convex  mathod  is  the  only  procedure  for  chancel  flood  routing 
available,  and  th:J  ptorage-'indicat ion  msthod  is  the  only  procedure;  for 
reservoir  flood  routing  availablac    Thoss  are  mathodo  ubich  ars  easily 
adaptable  to  Ey.nual  solution  and  perhaps  nere  chosan  for  thi&-  reason* 
Possibly,  other  methods  would  be  more  realistic  ard  as  easily  handled  by 
the  computer..    By  sinply  adding  other  routing  subroutines,  the  overall 
couiputcr  model  could  bo  wB.de  more  flexibJ.Gc 

Design-djvpa  Stormf^  ■  • 

An  example  of  design  peak  discharge  determination  based  on  frequency 
of  storm  was  given  in  Chapter  VII  ard  needs  no  further  amplification  hare* 
It  might  be  pointed  out,  bonever,  that  the  accuracy  with  which  the  i)Qd.\i 
discharge  of  a  hj'pothetical  stoj.n  can  be  predicted  can,  be  no  better  than 
the  accuracy  of  the  determinat ioxi  of  the  amount  of  rainfall  to  be  ex- 
pected in  any  given  location  for  the  storm  frequency/  chosen, 

Tba  hydrographs  which  were  detorminad  for  design  stormo  should  ha 
considered  only  for  their  quaj.itative  meaning  and  not  their  quantitative 
results c    The  reason  for  this ,  of  course,  is  that  the  model  of  the  water- 
shed  was  not  thoroughly  adjusted  in  the  sense  describod  in  Chapter  III. 
Furt?j.er  adjustments  in  the  model  or  watershed  param^^ters  could  now  be 
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■mde  quite  readily.,    With  a  DatisfactoriJ.y  adjvisted  model  available^  the 
hypothetical  stcriiiG  could  ba  routed  through, the  model  once  more,  and  de- 
sign hydrography  xjoulcl  he  tharefcy  obtainedo    Ac  furuhor  runoff  events 
occur  on  Back  Creek,  it  iray  bsconie  possible  to  determine  hou  the  runoff 
pattern  changes  v?ith  oizs  of  storiai  thif?  V;0ulcl  permit  iiiodtf ications  to 
be  Iliads  to  the  viatershed  model,  uhich  could  then  agaiD  be  sijhjceted  to 
the  hypothetical  stoms,  and  still  more  definitive  results  obtained. 
EE2?c^  ■^■'^~i£i^i!L^  State  High}-?ay  ^P^^^^i^t 

■  The  results  of  the  study  indicate  that  before  the  Konta;ea  State 
Highway  .Depart?iic;nt  could  reliably  use  the  SCS  fiat  hod  of  predicting  pro- 
bable peak  discharge  from  Montana  watersheds a  long  rarge  study  and 
collection  of  data  uould  first  be  neceosaryc    Befox'e  the  SCS  method  could 
be  applied,  liiaps  locating  soil  types  and  land  use  on  the  watershed,  or 
a  Eir;g3.e  i"ap  delineating  the  soil-cover  eo-uplex  numbers  (CjJ)  on  the  water- 
shed, would  have  to  be  made  available «     land  use  maps  couJd  be  drawn 
largely  i?ith  the  use  of  aerial  photographs  and  a  miniiiial  firsthand  knovj- 
ledge  of  the  location  involved^    Fydrologic  soil  types  can  presently  be 
deteriiiinad  in  many  areas  in  Montana  from  e?^istiiig  data  from  SCS  soil 
sui-veySo     Ho'ievcr,  large  sections  of  I^ontana  have  not  yet  been  surveycdc 
Therefore,  a  complete  soil  type  isap  for  the  state  could  not  now  be  com- 
piled» 

For  prediction  of  a  design  peak  discharge  fron  a  particular  watershed, 
additional  data  vould  have  to  be  collected^    Areas  of  subv'atersheds  and 
chamel  reach  lengths  could  be  i&easured  froru  aerial  photographs  if  they 


V 


were  available  for  the  uat-ershed  in  question*,    B'ield  Euiveys  v;ould  be 
necessary  to  determina  chax^nel  cross-sections «     Channel  bed  slojas  vjould 
also  have  to  ba  dctornilned  in  the  field  unless  suitable  topographic  mps 
were  available «    The  Manning  roughness  coefficient  could  be  determined  by 
inspecting  the  area  or  by  vie'n.rg  photographs  of  the  charnelo 

■  Kbntana  State  University  is  presently  installirg  a  third-generat ion 
computer;,  the  S13  Sigma  7,  vhicb  will  be  capable  of  running  the  current. 
(1967)  version  of  the  SC3-PFP  with  a  fey  modifications e  Unfortunately, 
the  Sigiiia  7  ^"2S  not  available  during  the  tiiue  of  this  study ^  aBl  it  '^as 
not  feasible  to  make  more  than  the  tvo  sets  of  runs  previously  described. 
To  effectively  make  use  of  such,  a  prcgram,  the  hydrologis g  should  have 
easy  access  to  the  cci'ipuber  because  of  the  ssveral  trials  uhich  seem  to 
be  necessary  in  adjustmenb  and  use  of  a  modele 

In  sur&iary,  the  SCS  method  requires  collection  of  a  reasonably  large 
amount  of  data ,     It  has  not  been  pass ible^  with  the  data  uhieh  were  avail- 
able for  this  study;  to  evaluate  the  SCS  mrethod  to  determine  v/hether  such 
amounts  of  data  collection  are  justified^. 


Cbap-fcer  IX 

coi^cnjsicis 

Ifce  investigation  reportad  bars  In  bar?  been  a  study  of  tlio  SCS 
method    of  predicting  peak  discbarge  flov  rates  £rori-  rain-caused  storas 
on  small  watersheds.    Data  for  one  actual  storm  which  occurred  on  Duck 
Creak  in  Prairie  and  KcCone  Count ie;?^  Montana together  vjith  eavaral 
bypotbet j-cal  stopins  for  the  same  watershed,  were  processed  hj  a  coiri- 
puterlsed  solution  of  the  SCS  method ,    The  resulting  synthet ic  hydro- 
graph  for  tb-^  actual  storm  was  compared  with  the  actual  hydrograpb  which 
was  recoru:^d  for  the  event*  ■ 

Findings  fro;;!  the  study  lead  to  the  following  conclusion's: 

1  "  The  SCS  method  is  a  systerBtic.  logically  organized  procedure 
which  generates  a  synthetic  bydrograpb  for  a  srall  watershed,,'  It 
utlj.izes  data  characterising  the  causative  rainstorm,  watershed  para- 
meters, and  accepted  flood  routing  tech^^iques.     The  iDetbod  has  been 
effectively  prcgraciiaied  by  the  SCS  for  coiuputer  solution-, 

2  -  Coi'iparison  of  the  synthetic  hydrograph  gensi^ateo  by  ti^e  coinpirber 
program  for  the  storm  on  June  l6,  l^o^?  on  Dusk  Creek  with  the  actual 
hydrograph  which  was  recorded  for  the  sane  event,  shows  considerable  dis- 
crepancy. This  is  thought  to  be  priiviarily  caused  by  using  iri^proper  value 
for  watershed  parameters  and  flood  routing  coefficients  as  input  to  the 
computer  program.  Alternative  factors  which  nay  explain  part  of  the  dis- 
crepancy iiiclude  the  possibilj;ty  that  the  actual  hydrograph  for  the  event 
fi:ay  not  be  entirely  accurate,  and  the  possibility  that  certain  tables  and 
graphs  utilized  by  the  SCS  not hod  nay  not  reflect  the  proper  values  wblch 
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should  be  used.  In  I-bntam. 

■  3  °  Suco^Gsful  use  of  the  SC3  iri^thod  rcQuir^so  carefiii  defi2>ition  of 
uatorsbed  Z3.T8.u3-te:cB  and  flocd  roiri;i:-g  coefficients^    Vfhen  a  satxsfactoi'y 
model  is  finally  obuaiDed,  ths  netbod  represents  a  very  effec«5.v&,  rapid 
method  for  deter;nird.Dg  paak  flow  rates  to  be  expscted  from  an  actual  or 
'  h;ypot  bet  leal  storm^  ;        ■        .  . 


Chapter  X 
RECOMMSNDATIOrB  AKO  SUlC-IARY 

Since  the  t:^c:U^l  of  Diicl^  Croe:k  watershed  has  mrj  hoan  constz-actcd 
anci  partially  ad,-]u3tt?d,  further  adjusta^ant  should  be  r.mde  as  Kore  run- 
off events  occur  on^  Duck  Creeks    Also_,  wh^n  additlOijRl  discharge  Reasure-- 
ments  are  iBds  by  the  USGS^  the  reliabilj.ty  of  the  prescmt  rating  curve 
at  the  I^ontana  State  University  rscordor  on  Dack  Creek  should  be  checked ^ 

Models  of  other  Montana  :^'at;;;rsheds  chouri  be  con-^trrcttd  to  ascer- 
taizi  the  variability  in  i-u'^off  characteristics  aruong  them*    Duck  Croak 
isay  not  be  representative  of  eastern  Montana  vatershGds  for  soius  unfore- 
seen reasons .    Thus ,  v/ith  other  inodels  available,  the  errors  in  design 
storm  peak  discharge  predictions  would  be  diuiinlshed^ 

Another  area  of  potential  study  vould  be  the  testing  of  th?  various 
facets  of  the  SC3  msthodc    Since  the  present  program  is  based  largely  on 
the  method  used  before  the  advent  of  the  digital  cov.iiputer  program,  there 
are  cciriponents  of  the  nssthod  that  could  be  iisre  sophisticated  without 
making  the  problem  solution  Hore  difficult  ^,     It  uould  appear  that  the  syn^ 
thetic  hydrograph  shape  could  be  iiade  Dore  flexibJ.e  to  account  for  dif- 
ferences in  area,  shape  of  vatershed,  cliiratic  condition  and  lGcaticn„ 
Calculation  of  the  thsoretlcal  tise  of  concentration  and  lag  time  needs 
f mother  study.     The  iiuportanca  of  syjruhesij^ing  the  tiining  of  runoff  has 
been  sVioun  by  the  cxav-ple  used  herein. 

In  order  to  PS3.ke  fu.rtber  study  both  easier  and  more  profitable;  the 
SC3-FPP  should  be  altered  slightly  for  use  on  the  rev;  digital  co-rputer 
presently  being  instslled  at  Montana  State  University, 
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Suiraiiary 

Prediction  of  the  peak  discharge  of  a  rain-causad  runoff  event  from 
a  watershed  by  the  n'ethod  developed  by  the  SCS  has  been  discussed.  The 
developmsnt  of  this  iaatbod  vjas  outlined  in  Chapter  III<, 

I^fere  recently  the  i^iothod  hss  been  prograimnGd  for  use  with  an  elec- 
tronic eosiputero    A  description  of  this  prcgrara  and"  an  explanation  of  its 
operation  Ir;:  fonnd  in  Ch-'ipter  IV*  ' 

The  ability  of  tha  SCS  method  to  accni'ately  predict  peak  discharge 
was  partially  tested  with  use  of  the  computer.    The  nodal  of  the  Duck 
Creek  vmtershad  iiaz  constructed  and  \ia.s  subjected  to  the  storr^  of  June  l6^ 
1965^    Although  the  test  run  is  not  Kufiicient  to  evaluate  the  predict in.g 
ability  of  the  prO{.^ram,  the  model  is  nov?  available  for  Duck  Creek  and  can- 
be  further  adjueted  by  testing  other  storms. 
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APPENDIX 


Direct  Runoff  (Q)  in  Inches 
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APPENDIX  B 


Hydrologlc  Soil  l^P^ 

A.  (Low  runoff  potential).    Soils  having  high  infiltration  rates 
even  when  thoroughly  wetted  and  consisting  chiefly  of  deep,  well 
to  excessively  drained  sands  or  gravels.    These  soils  have  a  high 
rate  of  water  transmission. 

B.  Soils  having  moderate  infiltration  rates  when  thoroughly  wetted  and 
consisting  chiefly  of  moderately  deep  to  deep,  moderately  well  to 
v?ell  drained  soils  with  moderately  fine  to  moderately  coarse  tex- 
tures.   These  soils  have  a  moderate  rate  of  water  transmission, 

C.  Soils  having  slow  infiltration  rates  when  thoroughly  wetted  and  con 
Bisting  chiefly  of  soils  with  a  layer  that  impedes  downward  move- 
ment of  water,  or  soils  with  moderately  fine  to  fine  texture.  Thes 
soils  have  a  slow  rate  of  water  transmission. 

D.  (High  runoff  potential).    Soils  having  very  slow  infiltration  rates 
when  thoroughly  wetted  and  consisting  chiefly  of  clay  soils  with  a 
high  swelling  potential,  soils  with  a  permanent  high  water  table, 
soils  with  a  claypan  or  clay  layer  at  or  near  the  surface,  and 
shallow  soils  over  nearly  impervious  material.     These  soils  have  a 
very  slow  rate  of  water  transmission. 


APPENDIX  C 

Runoff  curve  numbers  for  hydrologic  soil-cover  complexes 


(Antecedent  moisture  condition  II,  and  I^  =  0,2  S) 
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APPENDIX  D 

MANNING  Q  CALCULATION   FROM   X-SFCTION  NOTES 

>>. 

C  A  =  AREA  BETWEEN   TWO   X-FECTIorj  POHlTS 

C  AX3EC  =   TOTAL  AREA  OF  E.-ITIRE  X-SECTIO.N 

C  ATOTL  -~   TOTAL  AREA  OF  PORTION  OF  X-SECTION 

FOR  EXAMPLE,       AREA  OF  LEFT   FLOOD,  CHANNEL, 

OR  RIGHT   FLOOD  PLAIN 

C  COEFN  =  MANNING     W  CCEFFICIE^IT 

C  DELWS  =  DISTANCE   Y  TO  WATER  SURFACE   FROM  r!EXT  POINT 

C  (X,Y)    ABOVE  WATER  SURFACE 

C  IHWSEL  =  HIGHEST  WATER  SURFACE   TO  BE    I  NVE.^T  I  GAT  ED 

C  LCWSEL  ^   LOWEST   '^'ATER   SURFACE   TO  BE  INVESTIGATED 

C  NCWSEL  --   NO.    OF  WATER   SURFACE  ELEVATIONS  TESTED 

C  NOXSA  ^   iiC.    OF   X-SECTION  AREAS 

C  NOXSP  =   NO.    OF  X-SECTIOfl   POINTS,    10    IS   THE  MAXIMUM 

C  NSU3XS  =   NEW  SUB  X-SECTIOr; 

C  IF  ^   0     NEW  SUB-SECT  RW]   TO  BE  READ 

C  IF  ^    1     NEW  X-SECTION   TO  BE  READ 

C  .       P  =  PERi:'ETER 

C  FTOTL  ^   TOTAL   PERI'-IETER    IN   SUB  X-SECTION 

C  PX5EC  PERIMETER'  OF   WHOLE  X-SECTION 

C  SAVXL  ^   LOCATION  TO  SAVE  VALUE  OF   X  ON   LEFT   SIDE  OF  WS 

C  Ii.!TERSECT  !0:.l   TO  OF  HSED  WITH  NE'-'  ^JS  ELEVATIONS 

C  SAVXR  ^   SAME  AS   SAVXL   EXCEPT   FOR  RIGHT  SIDE 

C  SLOPE  =   CHAr.riEL  SLOPE 

C  0  =   TOTAL  DISCHARGE   IN  SUB  X-SECTIOi>; 

C  QXSEC  =   DISCHARGE   i^ROM  WHOLE  X-SECTION 

C  VXSEC  -   AVERAGE  VELOCITY   FOR  WHOLE  X-SECTIO,N 

C  WSELEV  =   WATER   SURI^ACE  ELEVATION 

C  X(I)  =   DISTANCE  FROM  HORZ  CONTROL  ■ 

C  Y{I)  -   ELEV.   ABOVE  DAT!JM    '  ■      '  ■ 
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C  DISCHARGE  Q  CAM  EE   CALCULATED   FROM  MAMniflGS   EQUATION  FOR 

C  IRREGULAR  f>lATURAL  STREAM  CRO  5S- SEC  T  I  Oi  OS  .      THE  WHOLE  X-SEC 

C  MAY  BE  SUBDIVIDED    UNTO  SUB-SECT  I  QMS   AT   SHARP  GEOMETRICAL 

C  BREAKS.      THERE  ARE   FOUR  DATA  CARDS  .MECESSARY  FOR   EACH  SUB- 

C  SECT  I  on.        THERE  FORMAT    IS  AS  FOLLOWS. 
C  ■  '      ■  , 

C  CARD  COLUMMS 


c 

1 

30 

PRINT 

Air^  TITLE  ShIOWING  WATERS 

HED 

c 

LOCAT I  or;,    SECTION  NUMBER,  ETC. 

c 

-  2 

1 

3 

f^OXSP 

RIGHT  JUSTIFIED 

c 

6 

10 

LOWSEL 

RIGHT  JUSTIFIED 

c 

1 1 

15 

I HWSEL 

RIGHT  JUSTIFIED 

r 

\- 

1  6 

20 

NSUBXS 

RIGHT  JUSTIFIED 

c 

21 

3C 

COEFN 

I      FORMAT      X  X  X  X  .  X  X  X  X  X 

c 

31 

SLOPE 

I  \\    F 0 R AT      X  X X  X  .  X  X  X X  X 

c 

3; 

1 

G 

X  DIST 

TC   15T   POiriT   FROM  REFER 

ENCE 

c 

'•MTH  FORMAT  XXXXX.XX 

9 

\L 

X  DIST 

TO   2f!D   PT   FROi'l  REFEREf^C 

E  PT 

c 

17 

24 

.3RD 

c 

ET 

r 

c  . 

/, 

1 

8 

ELEVAT 

ION  Y  CORRESPONDING  TO  X 

OF 

ETC. 

C  AFTER  ALL   SUB-SEC  T I  OrJS  HAvE  RUN  FOR  A  GIVEN  X-SECTION 

C  ANOTHER  CARD  MUST  BE   READ   TO    LNITIALIXE   THE  SU'-IMARY  PRINT 

C  OPERATION.      THE  CARD  MUST   HAVE   TfiE  LOWEST   VALUE  OF  LOWSEL 

C  FOR  THE   X-SECTION   LN  COLUMN    1-5',    RIGHT  JUSTIFIED. 
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D  I  MEMS  I  0[-J   T  I  TLE  (  4  '  )  ,  X  (  1  J  )  ,  Y  (  1  0  )  ?  [^XD  Y  (  9  )  ,  A  (  9  )  ?  P  (  9  )  *  D  (  1  0  ) 
1  VXSEC  (  6  -  )  5  AXSrC  (  5n'^  )  ,PXSlC  (  600  )  ,QXSEC  (  &'"-0  ) 
PRINT  104 

104  FOPMAT  (  oOfjO  MANNITIG  0  CALCULATION  FROM 


1 X-GLCT I 


r 


C  iNOTE   TMA^    THE   FOLLOWING  DO   LOOP    IS   GOOD   OMLY   -OP  ELEVA- 

C  TlOliS  FPO'i  2  300   TO   3100   FEET     CHAriCE  ELEV.    HI  STATEi'^EMT 

C  19   TO  FIT   WATERSHED, DIMEMS I  C:j  STATEMENT   MUST  ALSO  BE 

C  CHAi'jGED   FOR  VXSEC,    AXSEC,    PXSEC,    AMD  QXSEC.  CONSTANTS 

C  IN   STATEMENTS  2^-,    105,   A[>]D    107  MUST  CONFORM 

19  DO   2  3   k;<^2  30  1  ,  3  1  00 

2  4  K-KK~25JJ 

AXSEC  ( ;<  )  -  0.0 
PXSEC  CO  ^-0  .  0 
.23   QXSEC  CO  =0.0 

1  READ   89,  ( T  ITLE(  I  )  , 1  =  1  ,40) 

3  9   FORMAT ( 40A2 ) 

PRINT   90  ,  ( T  I  TLE (  I  )  ,  I  =  i  ,40) 
^     90   FORMAT ( IHl ,40A2/// ) 

READ  91  ,NOXSP  ,LOWSEL  ,  I  HWSEL  ,  NSLiBXS  ,COEFN  ,SLOde 

91  FORMAT(4I5?2F10.3) 
IF(iNOXSr^)  16,16,2 

2  READ  92  ,  { X {  I  )  , 1= 1 ,NOXSP)  -  '  ' 
READ  92  ,  (  Y(  I  )  ,  1  =  1  ,,NOXSP)                   "     ,  . 

92  FOR.": AT  (  iOF3. 2  )  ' 
PRINT    93 ,  ( X(  I  )  ,    1  =  1  ,NCXSP) 

93  F0RMAT(4H   X=    l-^FS'.O/)  _  ' 
PRINT    1^3 ,  ( Y(  I )  ,  1  =  1  ,NOXSP ) 

103   F0RMAT(4FI   Y=  10F0,2//) 

N0XSA=N0XSP-1 

PRIiiT   95,  COEFN 
93   FORMAT(9H0         N     =F10.5/)  :     "      ■  - 

PRINT   97,  SLOPE 
97  FORMAT (    9H  SLOPE  =F10.3///) 

PRINT  101 

101   FORMAT (62H  WATER  SURFACE       DISCHARGE  Q     END  AREA 

1  VELOCITY       PERIMETER       /61H         ELEVATIO.N  (CFS) 

2  ( SO   ^T )  { EPS )  { FT  )  /  /  )  _ 
DO    1  3   J   =LO';;SEL  ,  I  H^-ZSEL 


,1 


TOTL=0. 


PTOTL=0.0 
DO  3  1=1, 
AJ  =  J 

D(  I  ) =AJ-Y 
IF    (0(1))  , 

4  D  (  I  )  =  0 

3  CONTINl 
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DO   1  o   I  -  1  ?  I  ^0  X S  A 
I  I  ^  1  + 1 

I  F  (  Y  (  I  I  )  - (  I  )  )    5  ,  1  7  ,  5 

1 7  DXDY ( I ) -0 
GO   TO  18 

5  DXDY (  I  ) - ( X (  I  I  ) -  X {  I  )  ) / ( Y {  I  I  ) - Y (  I  )  ) 

18  CCMTiriUE 

C  CALC'JLATl  X  AND   Y  C00RDIMAT[:5  WrIERE  INTERSECTS 

C  X-SECTICM  on  LEFT   SIDE   ^ACIi'G  DCV/NST  RE  AM ,        THE  IF(DXD^ 

C  BYPASSES   THIS   SIDE    IF   X-SECTICM    IS  OF  RIGHT   FLOOD  SEC- 

C  TlOr:  ONLY 

K-1 

IF{DXDY(!<)  )    6,  SI  ,31 
51   SAVXL  =  X(;<) 

GO   TO  9  - 

6  K=K+1 

IF    {D(K))  6,6,7 
-      7   KK  =  K-1  .  . 

SAVXL=X ( KK )  .    •  ■ 

DELV.'S=^Y(  !<K)  -AJ 

IF    (DELWS)    6,9,3  '  ■ 

'6   X(  KK  )  -X(  KX)  -DELW9^DXDY  (  ) 
C  CALCULATE  X  AND  Y  COORDINATES  '-/HERE  N'.S.  INTERESTS 

C  X-SECTIO:^!  01'!   RIGHT  SIDE   FACIf^iG  DOWNSTREAM 

9  L^riCXSP 

LLL-L-1  "  _  ■ 

I F (DXD Y ( LLL) )  32,82,10 
82  SAVXR=X(L) 
LL-L 
L^LLL 

GO   TO   13  ■  .    ■"   •  ■  ' 

10  L=L"1 

I  F    (  D  (  L  )  )    1 0  ,  1  3  ,  1 1 

11  LL-L+I  - 
SAVXR-X ( LL ) 

DEL'/S^Y(  LL ) --AJ  ■     .  • 

IF    (DELV.'S)  13,13,12 

12  X  (  LL  )  =X  (  LL  )  -DELV./5-"DXDY(  L  ) 

C  CALCULATE  AREA  ,  PERM  I  ME  TER  ,    AfJD  MAr^NHlG  0 

13  DO   14    I=:'K  ,L  ■■-  ■ 

II  -  I  +  1 

A  (  I  }    (  AD  S  F  {  X  f  I  I  }  -  X  (  I  )  )  )  -  {  (  D  {  I  I  )  ^    {  I)  )  /  2  .  ) 

P  {  I  )  - SOR T  F  (  (  (  A3  S  F  (  X  (  I  P  -  X  (  I  )  )  )  - 2  .  )  +  {  (  AD  S  F  (  D  (  I  )  -  D  (  I  I  )  )  ) 
l-"2,)) 

ATOTL^ATOTL  +  A (  1)  . 

14  PTOTL  =  PTOTL^-P  (  I  )  ' 
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Q=  (  I  .49/CGEF"i  (  ATCTL-^^l  .  5  67  )  /  (  P  TOT  0  .  6  6  7  )  )  -  (  SLOPE 
1  -  -  0  .  5  ) 

AVEVEL=Q/ATOTL 

X ( LL ) -SAVXR  ,    ■  ■ 


X(  K 


1/  1  -  Q  A 


AVXL 


PRirj    lO:^,      J  .0  ,  UOTL  ,  AVEVEL  .PTCTL 
100   FORMAT  (  3  X  ,  I  5  f  8  X  ,  r- 1    .    ,  3  X  ,    3  .  1  ,  3  y  ,  ^6  .  2  »  5 X  ,  ^ *  1  ) 

PXSEC ( L ) -PXSEC ( L) +PT0TL 

AX5EC  (  L  }  ==AX5EC  (  L  ) +AT0TL 

GX5EC  (  L) --^GXSEC  (  L  ] +0 
15  VXSEC ( L ) =OXSEC ( L ) /AX5EC ( L ) 

IF(ri5U3X^:)2C,l  r2C^ 
2  0   PR  I  .XT   1  D5 


105   FOR'XAT  (38!' 


SUMMARY  FOR   Ei]TIRE  CROSS   SECTIOM  /// 


PRir;T  10  1 
READ  22 


LOWS EL 


22    FORr'AT  (15) 

DO   21   J   ^L^WSEL  ,  I  H';'^EL 
107  L-J-2500 
21    PRIMT   13-, J    ,GXSEC { L )  5 AXSEC ( L ;  rVXSEC ( L 

GO  TO  19 

p  i_  A  G  E  T  V;  0  EL  A  M  K  D  A  T  A  C  A  RDS   A  T   E  X  D  0  ^  DAT  A 
15  CALL  EXIT 


n  V'  r  rr  /"  '  I  ) 


i_iMI 
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